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Abstract. In this paper, we extend the multiscale approach developed in [Abel et. 
al.. Rep. Prog. Phys., submitted] by exploiting the scale separation between ions 
and the electrons. The gyrokinetic equation is expanded in powers of the electron 
to ion mass ratio, which provides a rigorous method for deriving the reduced electron 
model. We prove that ion-scale electromagnetic turbulence cannot change the magnetic 
topology, and argue that to lowest order the magnetic field lies on fluctuating flux 
surfaces. These flux surfaces are used to construct magnetic coordinates, and in 
these coordinates a closed system of equations for the electron response to ion-scale 
turbulence is derived. All fast electron timescales have been eliminated from these 
equations. Wc also use these magnetic surfaces to construct transport equations for 
electrons and for electron heat in terms of the reduced electron model. 
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1. Introduction 

The energy, particle, and momentum confinement of present-day fusion experiments and 
proposed future devices is limited by turbulent, rather than coUisional transport. The 
turbulence that causes this transport has an essentially multiscale character. It occurs 
on (perpendicular) spatial scales smaller than those associated with the equilibrium and 
on timescales shorter than those associated with transport of mean quantities but much 
longer than the gyroperiod. This is the basis for the multiscale gyrokinetic approach 
presented in [1, 2], characterised by the small parameter e = p/a, where p is the thermal 
gyroradius and a is a typical equilibrium length scale. In these treatments, all species 
are considered equal. However, there is a second scale separation that is important, the 
scale separation between the plasma ions and the electrons, characterised by the small 
parameter 6 = ^^/me/rrii, where mg and are the electron and ion masses respectively. 
For a deuterium plasma 6 ~ 1/60. 

Ion-scale fluctuations (e.g. Ion Temperature Gradient (ITG) [3, 4, 5], Parallel 
Velocity Gradient (PVG) [6, 7] or Trapped Electron Mode (TEM) [8] driven turbulence) 
have length scales perpendicular to the magnetic field comparable to the ion gyroradius 
and have frequencies comparable to the ion transit or bounce frequencies. In contrast, 
electron-scale fluctuations have perpendicular length scales comparable to the electron 
gyroradius, and frequencies comparable to the electron transit or bounce frequency. 
Unless they are suppressed (e.g. by strongly sheared flows), ion-scale fluctuations 
with their larger "mixing length" usually dominate the transport - and are therefore 
considered to be more important. In this paper we focus on the electron response to 
such fluctuations. 

The electrons are much lighter than the ions, hence ion-scale fluctuations are long- 
wavelength compared to the electron gyroradius and low-frequency when compared 
to the electron transit or bounce frequencies. It is clear that in simulating ion- 
scale turbulence one does not wish to resolve the short electron timescales. This has 
motivated the development of models for the electron response to ion-scale turbulence 
that somehow eliminate fast processes due to electrons [9, 10, 11, 12]. In this paper, 
we derive for the first time a set of equations governing the electron behaviour in the 
presence of ion-scale fluctuations. These equations eliminate the fast electron timescale 
via a rigorous expansion of the gyrokinetic equations in 6 <^ 1, subsidiary to the 
expansion in e ^ 1. 

The fast processes arise from the streaming of electrons along magnetic field lines. 
In electrostatic turbulence simulations, this streaming is along the unperturbed mean 
field and easily handled. The electrostatic approximation only holds at very low /3 (the 
ratio of thermal to magnetic pressure, and a measure of the importance of magnetic 
perturbations). However, any future reactor design will operate with as high a /3 as 
possible. Indeed, spherical tokamaks (such as MAST [ ]) already operate at finite (3. 
Thus, we require an electron response model that incorporates the fluctuations of the 
magnetic field in a fully general way. As we shall see, if we consider only ion-scale 
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Table 1. Typical length and time scales in selected fusion devices. 



fluctuations then the topology of the magnetic field is preserved by the fluctuations 
(see Section 3.2). However, to determine the electron behaviour we will need to make 
assumptions about the initial structure of the magnetic field and the effect of other scales. 
In particular, whether the field is regular, with nested flux surfaces, or stochastic. We 
show in Section 3.3, that if the magnetic field were significantly stochastic (i.e. if the 
ion heat transport due to the stochastic field were comparable to the transport due to 
the fluctuating E x B drifts) then the electron heat transport would be larger than 
the ion transport by a factor of ^ 1. Such poor electron heat confinement is not 
observed in experiments, indeed it is common for the level of ion heat transport to 
exceed the level of electron heat transport in conventional tokamaks. Thus, we assume 
that any stochastic component of the magnetic field is small enough to be ignored (see 
Section 3.3) and therefore that the field lines lie on fluctuating magnetic surfaces. 

The structure of the rest of the paper is as follows. In Section 2.1, we introduce the 
electron gyrokinetic equation and attendant notation. In Section 2.2, we develop the 
formal expansion in 5 ^ 1 and order all physical quantities with respect to this small 
paramater. The key result of this paper - the gyrokinetic system with reduced electron 
equations - is derived via a systematic expansion in 6. This derivation is performed in 
Sections 3 and 4 with the results concisely summarized in Section 5. 

In Section 3.2, we use the expansion in 5 ^ 1 to prove that long-wavelength low- 
frequency fluctuations preserve the topology of magnetic field lines, generalising the 
results of [14, 15, 16, 17]. With the assumption of an initially regular magnetic field, 
this result enables us to introduce a coordinate system aligned to the exact magnetic 
field in Section 3. It is this set of coordinates that allows us, in Sections 4.2 and 4.3, to 
average over the fast electron motion along exact magnetic field lines and so close our 
equations for the electron distribution function. These equations are particularly simple 
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in the limit of strong collisions (z/g ^ u) which we derive in Section 6. 

The second major result of this paper is a simplified system of transport equations 
for electrons. This is presented in Section 7 and derived in Appendix H. These equations 
emphasise the particularly simple form that the transport fluxes take in the low-mass- 
ratio limit. We also present the transport equations for the coUisional electron model of 
Section 6, which take on an even simpler form. 

We draw the results of the paper together in Section 8. We make a detailed 
comparison of our equations to previously-derived models in Section 8.1, emphasising 
both the differences and the similarities. Finally, we conclude in Section 8.2. 



2. Electron Gyrokinetics and the Mass- Ratio Expansion 

We study the response of electrons to ion-scale fluctuations within the framework of 
multiscale gyrokinetics [1]. First, we will introduce this framework in Section 2.1 along 
with all attendant notation. Then, in Section 2.2, we formally order all quantities in 
the small parameter 6 = ^/m^Jrrii. 



2.1. Gyrokinetics in a Rotating Tokamak Plasma 

In this section, we introduce the assumptions and resulting equations of this formalism. 
The detailed derivation of these results can be found in [1]. 

First, we split all physical quantities into mean and fluctuating parts: 

B = B + SB, B = (b) , (1) 



e « 1, (4) 



turb 

E = E + 5E, E = (e) , (2) 

\ /turb 

fs=Fs + 5fs, ^F, = (3) 

where B is the magnetic field, E the electric field, fs the distribution function of species 
s and is some average over the fluctuations. 

We assume that the fluctuations obey the standard gyrokinetic orderingi : 
u \6E\ \6B\ 6fs Vs ps 

VLs k± \E\ \B\ Fs fi, a 
where u is the typical fluctuation frequency in the frame rotating with the plasma, 
Qg = ZscB/msC is the cyclotron frequency of species s with charge ZgC and mass 
in the mean magnetic field B = \B\, k\i and k± are typical wavenumbers parallel and 
perpendicular to the mean magnetic field, is a typical collision frequency, ps = VthJ^s 
is the thermal Larmor radius, Vth^ = ^/'^Ts/rris is the thermal speed of species s with 
mean temperature Tg and a is a typical equilibrium length scale. The mean quantities 
{B, E, Fg) are assumed to vary on the long length scale a and the long (transport) 
timescale te, the energy confinement time. This long timescale is ordered so that 

~ e'^s. (5) 

I In this paper, as in [1], the symbol ^ is used to mean "is the same order as" (with respect to the 
appropriate expansion) rather than the more usual "is asymptotically equivalent to" . 
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Typical values for the timescales and length scales in a tokamak are given in Table (1), 
which shows the very large scale separation, characterised by the small parameter e. This 
scale separation in time and space of the mean and fluctuating quantities motivates us 
to define the average (■)t^J.b in terms of temporal and spatial averages over intermediate 
time and length scales. This is done precisely in equations (14)-(18) of [1]. 

It is shown in Section 3.2 of [I] that if we introduce the cylindrical coordinate 
system [R, z, 0) and assume axisymmetry with respect to 0, the mean magnetic field 
can be written (to lowest order in e) as 

B = I{tlj,t)V(j) + VtljxV(f), (6) 

where the poloidal flux function ip is given in terms of the mean vector potential A by 

V^(i?,^,t) = /l^ = i?2A-V0, (7) 

and /(■?/', t) = R^B ■ V0 is the toroidal component of the mean magnetic field. We will 
also need the following alternate form of the magnetic field [18, 19]: 

B = Vil)x Va, (8) 

where a is a variable that labels different field lines within a flux surface. We cut the 
toroidal surface on the inboard side to force a to be single-valued. This cut forms an 
axisymmetric ribbon with one edge the magnetic axis and the other the plasma boundary 
or separatrix. Letting I be the distance along a field line, {ip, a, /) is a set of field-aligned 
coordinates for B. There are many different forms of field aligned coordinates and none 
of our results depend on a particular choice. Nevertheless, we will use ('0, a, /) for 
specificity. Axisymmetric functions are functions of ip and / but not a. Note that both 

z) and /(V') are mean quantities and therefore vary on the transport timescale. 

We assume that the mean electric field is large, E ~ {vthjc) B, which can be shown 
to imply that the plasma rotates toroidally with a velocity that is species-independent 
and whose angular velocity only depends on the flux label ip [-"] and the slow transport 
timescale: 

u = u{ip,t)R^V(f). (9) 

If we express the mean electric field in terms of potentials using Gaussian units and 
Coulomb gauge, 

E^-V^-lf, (10) 
then the scalar potential takes the form 

!f = <l>{'^,t) + !fo, V'o-e^, (11) 
and we have the following relationship between u{tp,t) and $(■?/', t): 

.W..t)=c-. (12) 
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With these results for the fields, our orderings imply that the distribution function 
fs of species s takes the following form: 

fs =F, + 6fs, (13) 
F, = Fos {HRs),es, t) + Fi, {R„ e,, /i„ a,t) + (e^/) , (14) 



(15) 



Let us explain this notation. We have followed [ ' ] and changed phase-space variables 
from {r,v) to {Rs,es, fis,^,(j)- 

^ - (16) 



Rs = r - —b X w, 

Es = ^rrisV^ + Zse{<!>{ij,t) + i^o) - Zse<^{ijl,t), 



2B 



(17) 
(18) 



where b = B/B is the unit vector along the mean magnetic field, the peculiar velocity 



W IS 



w 



V 



u = w\\b + wj_ (cos -i? 62 — sin-i^ei) 



(19) 



where Ci and 62 are two arbitrary orthogonal unit vectors perpendicular to the magnetic 
field, the direction of the parallel motion is a = tun/lwjjl, and we have defined 



moC 



rs (r, V, t) = ^(r, t) + -^R'v ■ V0 



(20) 



The quantity tp*{r,v,t) is equal to c/{Zse) times the canonical toroidal angular 
momentum and is thus conserved in an exactly axisymmetric system - i.e., in the absence 
of fiuctuations. In these variables, (14) splits the mean distribution function into the 
near Maxwellian equilibrium 

-1 3/2 



-£s/Ts(lp(Rs,t),t) 



(21) 



27rTsi^{Rs,t),t)_ 

and the neoclassical distribution function Fu, which describes large-scale 0(e) deviations 
from a Maxwellian. The function Ng is related to the mean density Ug by 



Us = Ng{tlj{r,t),t)exp 



(22) 



We now drop the explicit slow time dependence of all quantities - this is for notational 
clarity only and is not an assumption that the mean quantities are constant. 

The fiuctuating distribution function, given by (15), is composed of the Boltzmann 
response (the first term of (15)), where the fiuctuating potentials S(p and SA have been 
combined to form the electrostatic potential in the frame rotating with the plasma 

S^' = s<^--u- 6A, (23) 

c 
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and the gyrokinetic distribution function hg. Assuming that all mean quantities 
are known, then the gyrokinetic distribution function hs is given by the gyrokinetic 
equation [1, 2]: 

dh 



dFos ^ rrisFos 



dip 



B 



du 

dip 



(24) 



where the guiding-centre drift velocity is 



-uj\ip, t)RWR - 2w\\uj{ip, t)b xWz + —V(po 

rria 



(25) 



and where we have defined the gyrokinetic potential 

X = — -I' ■ 5 A = 6ip' — -w ■ SA, 
c c 

and the fluctuating velocity due to x§ 

c 



(26) 
(27) 



We have also introduced the gyroaverage at constant -R^, Eg and fig- 

{x{r,w,t))j^= di3x{r{Rs,eg,fis,^),w{Rs,es,fis,^,cr),t) . (28) 

The equation for hg is closed through Maxwell's equations for the fluctuating fields. 
The fluctutating potential Sip' obeys the quasineutrahty condition: 



(29) 



where the integral over velocities is performed at constant r and the gyroaverage at 
constant r, w^^, and w± is defined by 

Ztc 

The fluctuating magnetic field is determined from 6A\\ = b-6A and 6B\\ = b-6B, which 
obey the parallel component of Ampere's law 

Air 



ddhg{Rg{r, w\\,w^,§),eg{r, w\\,w±,{}), ^g{r, w±), a, t). 



(30) 



Vi 5 A, 



^ZgC jd^ww\\{hg)^, 



(31) 



and the perpendicular part of Ampere's law (see the discussion in Section 7.5 of [1]) 
. 5B\\B 



47r 



V_lV_l : ^ j d-^w{mgW±w±hg)^ = 0, 



(32) 



§ This can be shown to consist, physically, of the fluctuating E x B drift in the rotating frame, the 
motion of guiding centres along perturbed field lines and the fluctuating Vi? drift. 
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respectively. The evolution equation (24) and constituent relations (29), (31) and (32) 
form a closed system of equations for determining h^, 6ip', 6A\\, 6B^\ on the turbulent 
timescale. Typically, the system (24), (29), (31), and (32) is solved until statistical 
equilibrium is reached (a few nonlinear turnover times). The slowly- varying equilibrium 
is treated as constant during this evolution. The slow time dependence of Fqs is 
determined by transport equations for uj{ip,t), Ns{ipyt)y Ts{ip,t) and I{4',t) and the 
equilibrium condition determining ip[r,t). Such equations are given in [1]. 



2.2. The Mass-Ratio Expansion 

In order to express the scale separation of the electron and ion scales, we introduce the 
secondary small parameter 5 = ^J^mjrrii <^ 1. We will use a subscript i to denote 
any ion species, of which there may be many. We consider the expansion in 5 ^ 1 
to be a subsidiary expansion of the gyrokinetic system (equations (24), (29), (31), and 
(32))||. Thus, we assume that e = pj/a, i.e., that the fundamental gyrokinetic expansion 
parameter is that of the ions, and that 

e < 5 < 1. (33) 

All other dimensionless parameters are treated as finite - i.e., independent of 5. 
Therefore, we are assuming that Tj ~ Tg and (5 = Snp/B"^ ~ 1, where p is the plasma 
pressure. In typical fusion plasmas, /3 is often small and one might worry about this 
assumption. A more detailed analysis shows that our expansion in 6 requires that [21] 

/3 » ^ = 5^ (34) 

which is equivalent to requiring that the Alfven speed is smaller than the electron 
thermal speed. ^ In plasmas of interest, this is indeed satisfied and we are justified in 
treating (3 as finite. 

We assume that fiuctuating quantities vary on ion rather than electron scales. Thus, 

we order the timescales of fiuctuating quantities in the rotating frame and the size of 

the fiuctuations by 

d d c 

— + w ■ V ~ vth,b ■ ~ — I (V^v^' X b) ■ V_l| ~ 

^1 (VM|| X h) ■ V^l ~ (V(55|| X b) . Vxl ~ ^ (35) 

so that the typical fiuctuation timescale, the ion parallel streaming time, and the ion 
nonlinear timescale due to the fiuctuating fields are all comparable. We assume that 
the toroidal rotation velocity u is comparable to the ion thermal speed u ~ Vthi, which 
implies that the electric field is ordered as 

E ^ ^B. (36) 
c 

II This is in contrast to the approach taken by [12], see Section 8.1 for details. 

^ For cases where /3 is so low that this does not hold, an expansion of slab gyrokinetics for very low f3 
plasmas has been carried out in [22]. 
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The perpendicular length scales are of course ordered as 

|Vx|~fc±~p-^ (37) 

Clearly, the turbulent dynamics of the ions are not simplified by the subsidiary expansion 
in S and hence hi is determined by the gyrokinetic equation (24). 
The electron dynamics, however, are simplified because 

hv,^. ~ "^^^ V ^ ^^^^11 X b) ■ Vxl » ^ + ■ V, (38) 
so the parallel electron motion is much faster than the evolution of the fluctuations and 
also 

k±pe ~ 5 < 1, (39) 

so we can neglect finite-electron-gyroradius effects. We choose the electron collision rate 
Ue to be comparable to the fluctuation frequency, i.e., 

i^e ~ — ~ (40) 

a 

Note that although electron collisions are comparable to the fluctuation frequency, and, 
therefore, important in the electron fluctuation dynamics, this would conventionally be 
considered a low collisionality plasma since Ues = i^id/vthi ~ ^e^/vthe ~ 5 ^ 1- Under 
these assumptions, we expand he, Sip', 6A^^ and SB^\ as 

he = /^W + /^W + . . . , 5^' = + + . . . , 

6A\\ = 6A\\ + 6A\\ + • • • , and = 6B\\ + 6B\\ + ■ ■ ■ , (41) 

where h^e^ ~ S h'^e^ etc.^ In Sections 3 and 4, we systematically expand the gyrokinetic 
system of equations, (24), (29), (31) and (32) in powers of 6 to obtain h^^ and close the 
system at first order in 6. 

3. Zeroth Order: Particular Solutions, Flux Conservation and 
Field-Aligned Coordinates 

In this section, we apply the orderings of Section 2.2 to the gyrokinetic equation for 
electrons and examine the consequences to lowest order. 

In order to expand the gyrokinetic equation, in Appendix A, we express 
explicitly in terms of the fields and find, (A. 3): 

(X>„ = - ^W) - +0 L'f) . (42) 



"I'-i o(|f) 0(4) 

+ In fact, we will only require dtp' or 6B^\ to lowest order in S. Thus, we will drop the superscripts on 
these fields. 
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Substituting this result and the expansion (41) into the gyrokinetic equation (24), we 
find, to lowest order in 6: 



^" ' dRe B [ ^ dR, j ' dR,, 

ew\\Foe 



(43) 



where SA^^ is evaluated at Re and we have used u{Re) ■ {dw\\/ dR^^ = (axisymmetry). 
These equations are written in terms of the variables -Re, ^e, A^e and i). It will be 
convenient to instead convert back to using r as our spatial variable. Since (43) is only 
accurate up to corrections of order 0(5e^f2jFoe) and all functions in (43) are evaluated 
at -Re, we can use the fact that k_i_Pf, ~ 5 to replace -Re by r everywhere. The function 
Ee is now given by the simplified expression: 

Ee = ^rriewl + fieB - C^q + O {d^Te) . (44) 

Defining b = -B/|-B| to be the unit vector along the exact magnetic field, we have 
5B_ 

13 ' ' " B 

Thus, after a little algebra (using (45)), we can rewrite all but one term of (43) as 
derivatives along the exact field line: 



6 = 6+^ + 0(e2) = b - i X VMi|(°)) + 0{5t). (45) 



(46) 



e. 



where we have used 



5fe = ^Foe + + O (e(5Foe) 



(47) 



the Maxwellian form (21) of Foe, and divided through by Fq^. Equation (46) is an 
inhomogeneous ordinary differential equation along the exact field lines for 6fe. 

In the remainder of this section, we will find the general solution of (46). The 
particular solution to the inhomogenous problem is constructed in Section 3.1. In order 
to solve the homogenous part of (46), it is necessary to discover the structure of the 
exact magnetic field. Thus, in Section 3.2, we prove that magnetic fiux is conserved to 
lowest order in 6 and hence that the topology of the magnetic field is fixed. This result 
then allows us in Section 3.4 to introduce field-aligned coordinates in which it is simple 
to solve the homogenous part of (46). This general solution will provide us with an 
equation for 6A» and place constraints upon the form of 
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3.1. Particular Solutions of (46) 

We can now find particular solutions of (46). To do this, we write 6fe in the following 
form: 

6fe = A(r, Ee, He, (^)Foe + ge{r, Be, /"e, 0") + O (e(5Foe) , (48) 

where Qe satisfies the homogenous equation 

w\\b ■ Vge = 0. (49) 

This can clearly be done for any 6fe by suitable choice of A. Substituting (48) into (46) 
we obtain 

^- - f ) ^ ^" ■ ^) - (I - ^) 

To find the form of A, we take the derivative of this equation with respect to /ig and the 
second derivative with respect to Se to obtain 

— = and b ■ V — 

respectively. Clearly, any solution to these equations is either independent of /ig and 
a linear function of Ee or satisfies b ■ VA = 0. Any solution of the second kind can be 
absorbed into the function g^, leaving us with the general solution for A given by 

where we have written the linear function of in the natural way as a perturbed 
Maxwellian. However, it should be noted that 6ne{r) is not the perturbed electron 
density and STe{r) is not the perturbed electron temperature because can have both 
density and energy moments. Inserting (52) into (50) and taking a derivative with 
respect to Eg, we find the following equation for STe'. 

b-V (\nTe + ^] =0. (53) 



b- V— = and b- = 0, (51) 



Finally, substituting both (52) and (53) into (50), we obtain: 

-V(|-f .-.)^^(|-v)m.. ,54, 

Familiar physics is contained in these two equations. It will be shown in Section 3.5 
that the contribution to g^ from passing electrons is a function only of tl) and not of 
a or /, and, therefore, so are its density and temperature. Thus, for the purposes of 
interpreting (54) and (53), we can add the density and temperature of the passing part 
of g^ to brif. and bT^, respectively. Hence, (53) describes how (part of) the temperature 
of the electrons is equalised along field lines, due to rapid thermal conduction. The only 
part of the temperature that does not equalize along the field line, as we shall see when 
we solve (49), is the part of the temperature of the trapped particles contained in g^. 
This is to be expected; the trapped particles do not traverse the entire field fine and so 
cannot equalise their temperature along it. 
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Equation (54) describes how the pressure gradient of the electrons (again, not 
including the trapped particle pressure) along the field line is balanced by a parallel 
electric field. As we shall demonstrate, it is useful to think of (54) as an evolution 
equation for 6A\\ where Sip' and Sue are given. 

3.2. Flux Conservation 

To solve (49), (53) and (54) we will need to know the spatial structure of the total 
magnetic field. As (54) determines SA\\, it determines 6B_i_ and thence b. In this section, 
we will prove that (54) implies that magnetic fiux is conserved and consequently that 
the evolution of the magnetic perturbation cannot change the structure of the magnetic 
field lines. 

Magnetic fiux is said to be conserved if there exists an effective velocity field ■UefT 
such that closed curves moving with this velocity always enclose the same amount of 
magnetic flux. It is proved by Newcomb in [ : ] that if such a Ues exists, it also preserves 
magnetic fleld lines and thus the magnetic topology is flxed. Now, if the parallel electric 
fleld satisfles 

Eb = -bVt (55) 
where ^ is a single-valued scalar function, then the velocity fleld 

u,s = ^(e + V^) X B + U{r)b, (56) 



52 

(where B = \B\ and U{r) is an arbitrary single- valued function) satisfles E + 
X B/c = — V^. Thus Faraday's law becomes dB/dt = V x (ugg x B). The familiar 
proof of the frozen- in theorem (that is reproduced in almost all MHD textbooks) clearly 
applies if we recognise Ues as the flux- and fleld-line-preserving velocity. The condition 
(55) is a special case of the general conditions for the existence of frozen flux that were 
investigated in [23]. 

In Appendix B, we prove that we can write E ■ b in terms of potentials as 

Eb=-bV^o-b- VV --^■b~-(^ + u-v] 6A:\ + O fe^^^) . (57) 

c ot c \ot J V c / 

In Appendix C, we prove that the mean magnetic flux is conserved irrespective of the 
fluctuations, and so (see (G.7)) 

.b = b-V^^, 58 

at eue 

where riie is the 0(e) correction to the mean electron density (see discussion before 
(C.7)). Finally, using (54) and (58) in (57), we obtain 



Eb = -bV 



H In A^e + V^o 



erif. e 



(59) 
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which is indeed in the form (55). f Thus, we can define 

J-^flnA'.+ ^^V^o + f. (60) 



where ^ is an arbitrary single-valued function that satisfies 

b-V|=0, (61) 

and with this definition, (56) defines the field-line-preserving velocity Ues- In (56), we 
are free to choose U (r) and ^ subject to the constraint b ■ = 0. We postpone making 
choices for these parameters until Section 3.5. Since the field is frozen to Ues the ion- 
scale turbulence cannot alter the magnetic field structure on the turbulent timescale. 

3. 3. Structure of the Magnetic Field 

In the previous section we have demonstrated that ion-scale turbulence preserves field 
lines. The topology of the field is thus fixed. But what is this topology? Electron- 
scale fiuctuations can alter the topology and so can generate a stochastic field| from 
an initially regular one. Such fiuctuations may arise "locally" from electron-scale 
instabilities (e.g. Electron Temperature Gradient instabilities or Microtearing) , but 
electron-scale fiuctuations may also be driven "non-locally" by a cascade from ion-scale 
fiuctuations. In both cases, some tanghng of the magnetic field is possible, indeed 
probable. Thus, at some intuitive level, the magnetic field is expected to be stochastic. 
A stochastic field will be advected by Uefr and remain stochastic. 

If the magnetic field were stochastic and stationary then the heat diffusivity of the 
electrons would be given by the well-known Rechester-Rosenbluth formula [24]: 

Xn = Dstvth, ~ vthjc ^^^1^^ , (62) 

where D^t is a species-independent diffusion coefficient characterising the stochasticity 
of the field§ and Ic is the parallel correlation length of the stochastic component of the 

f Strictly speaking, we have proven that E ■ b = b ■ + £'||6 where ^ e^VthiB/c. Thus the 
electric field is precisely E = — itcff x B/c — + E^^b. Substituting this into Faraday's law, we find 
that 

^ = V X (u^ff X - cV X . 

Estimating the size of the left hand side and the second term on the right hand side, we see that the 
correction £'|| will affect neither the lowest-order mean field B nor the lowest-order fluctuation 5B. 
Thus, we can safely ignore it. 

I A stochastic field is one in which points on neighbouring field lines become exponentially separated 

as they travel along field lines [24]. 

§ The field-line diffusivity Dgt is defined as [24] 

Dst = lim {6r{l) - 6r{0)f/2l, 

where Sr(l) is the radial displacement of a point after following a field line for a distance I and the 
overbar denotes statistical averaging over field lines (the factor of 2 in the denominator is chosen to 
agree with the definition in [''-!]). 
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field 5-Bstoch- Let us assume the stochastic component of the field is comparable to 
the ion-scale field perturbations (i.e., 5i?stoch/-B ~ e) and the correlation length is the 
system size (i.e., Ic ~ a). Then, (62) predicts that the electron heat diffusivity in a 
tokamak should be larger than the ion heat diffusivity (both the electrostatic and 
electromagnetic parts of the ion heat diffusivity). This is not backed up by experimental 
observations; ion and electron heat diffusivities are observed to be of the same order 
of magnitude. Thus, we take this as experimental justification for assuming that the 
amount of magnetic field stochasticty is limited. || 

We therefore assume good flux surfaces in solving Eqs. (49), (53) and (54) - these 
surfaces will be approximations to the exact field that we will refer to as the "ion-scale 
field". The ion-scale turbulence will distort the ion-scale field but the surfaces will 
remain intact. A smaller stochastic field component with Ic ^ a and 5B stock/ B ~ e5^/^ 
would cause observable levels of transport and we cannot rule out a component of 
this size. Indeed, such a field would change the final electron equations in this paper. 
However, we will assume for simplicity that the stochastic component of the field is zero 
to order 0{e5) (i.e., 6Bstoch/B <^ e6) in this thesis. 

3.4- Field-Aligned Coordinates 

If we write the ion-scale magnetic field in the Clebsch form 

B = Vip X Va, (63) 

then ip is a single- valued function that labels the flux surfaces and a is a function that 
labels different field lines within a given flux surface. To ensure that a is single-valued, 
we make a cut in the poloidal domain on the inboard side of the tokamak. This cut 
makes a ribbon forming a toroidal loop with one edge of the ribbon being the magnetic 
axis and the other the plasma boundary or separatrix. We allow the cut to move with 
the velocity Ucs- This is legitimate as u^s is continuous across the cut. Since the cut 
is perturbed by the turbulence, it will not now be axisymmetric. The expression (63) 
for the magnetic field allows us (see, e.g., [1 ]) to introduce a system of field-aligned 
coordinates a, I) where / measures distance along a given field line. By definition, 
/ runs from one side of the cut in the poloidal plane to the other side of the cut. As 

II Wc can formulate this condition precisely as a condition on the field-line diffusivity (defined in the 
footnote on page 13) . The limit on stochasticity can be expressed as 

where we obtain this estimate by insisting that the change in a function g following the field line is 
dominated by explicit variation along a field line {k\\g) rather than the field line connecting neighbouring 
points in the perpendicular plane {Dstk\g). The existence of a stochastic field giving rise to a field-line 
diffusivity of size Dst ^ epe would not invalidate the existence of 'ip and a, as (for the purposes of 
studying ion-scale turbulence) we could define ^ surfaces by the mean radial location of a field line 
averaged along a magnetic correlation length Ic- This introduces a displacement < of the ■0 surface 
from the exact position of the field line, but for turbulence where fcj_pe ^ 1 this correction is negligible. 
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V 



Figure 1. The exact field B and its ip = const, surfaces (red lines, blue surface) 
plotted atop the mean field B and a, ip = const, surface (black lines, green surface). 



SB <^ B, the variables aligned to the exact field only deviate slightly from those aligned 
to the background field: 

\ip — ipl eip, |a — a| ~ ea, |/ — /| ~ el. (64) 

Importantly, this means that all axisymmetric mean quantities are, to lowest order in 
e, functions of ip and I but not of a. 

We need expressions for the spatial derivative operators b V and the Poisson bracket 
{a, b} = b-(Va x V6) in our new coordinates. By using the fact that b-Vip = b Va = 0, 
we discover that 



6 . Va = 6 ■ V/ 



~ da 



dl 



da 
dl 



(65) 



for any function a{r). Similarly, by using the fact that b ■ (Va x Vip) = B to lowest 
order, we see that the Poisson bracket takes the form 

{a, b} = b-{Vax V6) + 0(e|Va||V6|) = B f - + 0(e|Va||V6|). (66) 

\dadip dip da J 
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Since ■Ueff is a field-line-preserving velocity, we have that 
d 



, Weff ■ V ) ^ = 0, (67) 



and 

d 



, lieff ■ VJ a = 0. (68) 

As / ~ / and V ■ u^s ~ eu, u^s also approximately preserves the length of field lines and 
we have, to lowest order in e, 
d 



+ Weff ■ Vj / = 0. (69) 
It is of course convenient to choose the cut in the definition of a such that it is convected 

with Ueff. 

Finally, the time derivative at constant ip, a and / is related to the time derivative 
at constant r by 
d 



dt 



= 1^ + **eff ■ V. (70) 



3.5. Solution of the Lowest- Order Equations 

We now use the magnetic coordinates introduced in Section 3.4 to solve (49) and (53). 
We split velocity space {e^ and /ig) into two regions: the passing region where, for a 
given ip, the equation 



W\\ (/, Ip, Be, He) = \/ — [Se " fJ'eB + e(fo] = (71) 

V rile 

has no solutions for any / (i.e., Ee > fieB — ecpo for all /) and the trapped region in which 
there are two values of / at which w\\ = 0. The points where w\\ vanishes are called the 
bounce points of the electron's orbit. In the passing region, (49) implies that 

^ = 0, (72) 
dl 

so Qe is constant along field lines. As most magnetic surfaces are irrational, one field 
line spans the entire surface and so, for passing electrons, cannot depend on a: 

ge = gpe{'4',£e,fJ'e,(^,t)- (73) 

In the trapped region of velocity space, (49) implies that 

either w\i{l,ip,ee, fJ^e) = 0, or — ^ = 0. (74) 

dl 

Thus, Qe must be constant along magnetic field lines, but only between the bounce 
points. So ge can depend on the location of the bounce points, and, therefore, on a: 

ge = gte{'ijj,0^,Be,fJ'e,t), (75) 

where the dependence on a has been dropped as ge must be continuous as passes 
through zero and changes sign. For convenience, we define gpe to be zero in the trapped 
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region and gte to be zero in the passing region, so the complete solution is just the sum 
of the two functions: 

ge = gpeii', £e, /ie, CT, t) + gte{^ , «, Ee, fie, t) . (76) 

The solution of (53) in our new coordinates is 

Te{llj)+6Te = fe{4^), (77) 

for some function T^. As we can absorb any part of ST^ that is only a function of ip into 
ge, we can pick Te^ip) = Te{'ip) and so 

We can also now pick the arbitrary function ^ in the definition (60) of ^ to be any 
function of ip. In Appendix E, we find convenient forms for the arbitrary functions Cii^) 
and U{r) to obtain (see (E.6)) 

u,s = u + ^bxV (V - C) , (79) 



where we have defined 



T 

c = - 

e 



80) 



ne V J dip 
With this definition, we can express the general solution of (46) as 

6fe =^i^Oe + (V^ + gpe{'ip,£e,fie,(^) + gte{'4',a,Ee,fie) ^^^^ 

+ O {eSFoe) , 
and the evolution equation (54) for SA^^ becomes 

+ w ■ M|| = cb ■ V (C - V) + O (e(5^],M||) . (82) 

This equation allows us to solve for 6A\\ correctly to lowest order in 6. Equations for the 
evolution of (, gpe and gte at the next order in our expansion are obtained in Section 4. 

3. 6. Solution for 6B± 

At this juncture, we have enough information to determine SB_i_ in two ways. Firstly, 
we have the evolution equations (67) and (68) for ip and a, with u^s given by (79). 
This enables us to calculate ip and a correctly to first order in e. In Appendix D, we 
demonstrate that this gives us enough information to construct 6A\\ from the equations 

^^-^ a,.d (83) 

dip dl da dl ^ ' 

Secondly, we could solve for 5A\\ from the evolution equation (82) and then use (83) to 
determine ip and a. Either of these methods is valid, and provides a solution for 5A\\ 
correct to the requisite order. For the purposes of simulating our equations, a third 
option presents itself: determine 5A\\ from (82), ip and a from (67) and (68), and then 
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use (83) as consistency checks upon the solutions thus obtained. We choose the first of 
these options, i.e., to eliminate 6A\\ in favour of the fluctuating fields — ip and a — a. 

One point remains to be resolved. It would appear that knowing ip and a determines 
6B completely (via (63)). However, to determine 6B\\ from ip and a, one must know ip 
and a to O(e^) or, equivalently, the compressible part of ■Ueg to higher order (0(e^UthJ; 
see Appendix D), but we do not. Therefore, we determine 6B\^ from (32) which involves 
only lowest-order quantities - quantities that we do calculate. 



3. 7. Parallel Electron Currents 

As explained in the previous section, we have obtained a solution for SA\\. However, 6A\\ 
is usually determined from the parallel component of Ampere's Law, (31). Substituting 
our expansion for (41) into the parallel component of Ampere's law (31), we find 

V^M|| = ^(y'd3«;^ll(/^(°) + /.«>^) -^XJ^^^ jSww\\{K)^, (84) 

where a sum over s = i denotes a sum over all ion species. In this equation, the term 
due to is larger than the rest by 6^^. Thus, to lowest order, 

j£ww\\hf^ = jSww\\gpe = 0. (85) 

This is a constraint on in the lowest-order solution given by (81). Note that the 
trapped particles do not contribute to the current constraint. The parallel electron fiow 
(and current) is comparable to that of the ions and contained in h'^e^. Defining the 
parallel electron fiow 

5u||e = jd?ww\\h^P (86) 

reduces (84) to an equation for 5u\\e: 

= 1^ V^^Ay +Y,ZJ d'ww^^ {K),. (87) 
s=i 

Thus, instead of determining 6A\^ from the electron current. Ampere's law determines 
the electron fiow from 6A\^ [21]. As we are no longer using 5 Ay to describe SB± in our 
system of equations, we consider V^(5A|| to be a shorthand for a complicated expression 
involving derivatives of ip — ip and a — a via the relations (D.IO) or (D.ll). We require 
6u\\e to complete the solution in the next order, but we do not need any information 
about /il^"* other than 6u\\p. 



4. First Order: Dynamical Equations for Electrons 

In this section, we return to the gyrokinetic equation (24) and proceed to the next order 
in the mass-ratio expansion in order to find evolution equations for and (. 
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Subtracting (43) from (24) and using the solutions for h'^'' and SA\^ from Section 3.5, 
we obtain the first-order part of the gyrokinetic equation: 



iMmewfiFQe duj , . ( e d^^l^ 1 r n^ 



where we have used (45), rewritten terms involving (V^)^ as Poisson brackets using 
(27) and (66), and, as in Section 3, replaced Re with r as our spatial variable. Thus, 
the collision operator in (88) is evaluated at Re = r. The drift velocity in (88) now 
contains only the leading-order part: 



Voe = --^^^ ymewf^b ■ Vb + -mewiVlnB - eV(poj + 0{e5vth,). (89) 
We now rewrite (46) in terms of 6fe and (: 



ImeWuFoe duj , X f e 95 A^}"* if , 

where we have used (47), (80) and the fact that 

Iteff- V = lt- V + |{V-C,-}, (91) 

which follows directly from (79) and (66). 

Finally, we substitute the form of 5fe from (81), to obtain 



dt 



Foe d_ 
Te dt 



_ (eC - fleSB^\) - Vj,e ■ V(^ - ^)^ + VDe " V (V - C) ^ 

'4>,a,l T e 



c d , , dFoe iWmewlFoedu , 



(92) 



+ w\\F( 



ll-roe 



e d5Af 
cTe dt 



+ l«.luF„.})+|«.Af} 
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where we have gathered all terms involving ge together on the left-hand side and all 
terms involving SA\^^'^ on the right. Notationally, (92) can be confusing as we have made 
some abbreviations for the sake of conciseness. Firstly, all spatial derivatives are taken 
at constant t, and /Ze- Thus, for equations that are solved in the moving coordinate 
system 



dip 



+ Va— 



+ vA 

dl 



(93) 



As this notation is cumbersome, we will continue to use V where appropriate (i.e., where 
it would expand to more than one derivative). Secondly, all the Poisson brackets in (92) 
are shorthand for the expression (66) in terms of ip and a derivatives. Finally, this 
equation is, intrinsically, an equation in the frame moving with the magnetic field, and 
all functions should be interpreted as functions of t, ip, <5, and / rather than t and r. 

Equation (92) describes how ge, and thus 6fe, evolves dynamically in response to 
given perturbed fields. As -UeS is also considered known, (67)-(69) can be solved to give 
the complete transformation from the fixed to the moving coordinate system. Thus, we 
consider all the quantities in (92) to be known in both the fixed coordinates (which will 
be needed when we come to the field equations) {ip,a,l), and the moving coordinates, 
{'ijj,a,l). However, in its current form, (92) is not closed, as we know neither h^e^ nor 
(JAii^-*. In the following sections we derive evolution equations for (, gte and gpe as 
solutbility constraints of (92). 



4.I. Fluctuating Continuity Equation 

The one piece of information that is known about /li^'* is its parallel velocity moment 
6u\u, determined by (87). Thus, we integrate (92) over all velocities and find 



dt 



T 



ne + / d^wge 



6Bii 



B 



-Up 



dl B 



c 



c refill 

e 



B 



d^Wflege 



^^°^ + (Vz, + v,)-v(V-C) 



d 



+ cne- [C-SB,^^ 



dijj 

T, \ fd\nN, 



(94) 



d\n.Tp 



dip 



dip 



d5B\\ n^TedhiTe 
da cB dip 



B diPdl 

where we have used Bh = B + 0{eB), the fact that h^e^ carries no current, and 

, ,o f BdSediied^ 

' d w ■- ^ ' 

We have also defined 



E 



V, 



D 



eB 



b X (b-Vb + VlnB) 



(95) 



(96) 
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and 



c 



b X Vv^o 



(97) 



B 



so that 



d'^wVBeFoe 



( 



Vd + K ) Ue. 



(98) 



We have also used V ■ -B = and (D.ll) to ehminate b and 6A\^ from our equation. 
It is important to note that this equation is written in the moving coordinate system. 



tjj, a and I. In this context, ip — ip should be viewed as just another fluctuating field. 
Therefore, all gradients should be thought of in terms of derivatives with respect to ip 
and a, for the Poisson bracket this is done via (66). 

That (94) is essentially the fluctuating continuity equation is not immediatly 
obvious. Let us attempt to make this more transparent. Clearly, the first two 
terms under the time derivative in (94) are the perturbed electron density due to the 
fluctuations; the time derivative of SB\\ adds to this the fluctuation in the electron 
density due to the compression of the flux surfaces.^ This total density is changed 
by various flows - the compression in the parallel flow 6ui\e, the flow of electrons past 
the flux surface given by {c/B)b x VC, the grad-i? drift due to SB\\, and the magnetic 
drifts. Indeed, even the final term in (94) can be interpreted as a compression - it is 
the compression of the electron fluid due to the perturbed flux surfaces "sticking out" 
into regions with different toroidal angular velocities. 

Equation (94) provides us with an evolution equation for ( if we know and the 
perturbed fields. The evolution equations for Qe given ( are found in the next two 
sections. 

4-2. Trapped Electrons and the Bounce Average 

In order to derive an equation for in the trapped region of velocity space, we need 
to define an averaging operator that eliminates the fast electron bounce motion (i.e., 
annhilates luyb ■ V and thereby eliminates /li^^). The appropriate average is the bounce 
average along the perturbed field i.e., at fixed ip and a. For an arbitrary function 
a{%lj, a, /, Be, /ie, cr) we define the bounce average to be 



where the integration is carried out from one bounce point where w\\ = 0, along the 
magnetic field, to the second bounce point where w\\ passes through zero and changes 

% In the derivation of the transport equations we show that 



It relates ({ip,a,l,t) to ge and the other fluctuating fields, also given as functions of 




(99) 




so d6B^\/dt is equal to the rate of perpendicular compression of the flux surfaces. 
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sign, and then back to the first bounce point. This integral is carried out at constant 
ip, a, Ee, fie and ^. Thus, equivalently, the bounce average is 

h ~ I h ~ 

(«)|| = [ [«(^ = 1) + = -1)] / 2 /" Ai, (100) 



'I 



\w 



III /J H\ 



where li and I2 are the bounce points of the particle orbit. 

Clearly, the bounce average commutes with the time derivative at fixed moving 
coordinates ip, a and /. In Appendix F, we prove that 

W||6- V/i(^)^^^ = (101) 
and that, to lowest order in e, 

yDe-V^)^^=0. (102) 

Using these results, we take (92) in the trapped region of velocity space and perform 
the bounce average to find 

d 



9te + (Voe ■ Vf7te)|| + ^ {e (C)|| - /^e ) ^ , } " ( ( C [(^te] | 



dt .-f'' ' '^^^'11 ' eB 

_ d_ 
~ Ft 



((C)||-/ie(55||>||)y^-(VDe-V(^-^ 



II dip 
dip 



(103) 



+ (Vne ■ V (V - C))|l ^ + ^ ^ (Oil - (5i?||>||) 

where the last line in (92) has vanished under the bounce average because it is 
antisymmetric in a (the sign of the parallel velocity) . This equation is a closed equation 
for Qte, once the fields and C, are given. Note that all quantities must be expressed in 
the moving coordinates, in order to perform the bounce average. Again, as with (94), 
all gradients should be interpreted in terms of derivatives with respect to the moving 
coordinates as should the Poisson bracket (via (66)). 

4-3. Passing Electrons and the Flux Surface Average 

In the passing region of velocity space, the fast streaming of electrons along perturbed 
field lines causes a passing electron to sample an entire perturbed fiux surface [ip = const, 
surface). Thus, the averaging procedure we will need is the average over a perturbed 
fiux surface. For any function a, this is defined by 



i^iy d^ra{r,ee,fie,^,(^) j J < 



(a(r, ^e, /^e, cr))^ = ^im^ <( j d-^ra{r,ee,iJ,e,i^,(T) / I d-^r } , (104) 



The Electron Response 



23 



where the region of integration A is the volume between the surface labelled by i/j 
and that labelled hj ip + Aip, and the integral is taken at constant e^, fie and i). In 
Appendix G, we prove that 



0, 



(105) 



that the flux surface average commutes with the time derivative following the magnetic 
field, that 



B 



\W\\ 



0, 



'4> 



and that, for any function a{r), 

' b X Va 



\Wl\ 



0. 



We also prove that, for any fluctuating function A, 



B 



\w 



Voe • VA 



0. 



(106) 



(107) 



(108) 



Taking (92) in the passing region of velocity space, multiplying by B/\w\i\, and 
flux-surface averaging (to eliminate w\ib ■ V/ii^''), we obtain: 




me|w;|||-^(^ - ^) 
ol 



^ Te dip 



Foe 



(109) 



InFf 




where the flnal term in (92) vanishes because in the passing region is only a function 
of ip and we can use (107). We still have to eliminate the flnal line of (109). We 
note that we only require the part of Qpe that is even in a - we require its density in 
the quasineutrality condition (29) and its pressure in perpendicular Ampere's law (32). 
Thus, we average (109) over the two signs of a and henceforth Qpf, denotes only the even 
part. Hence we obtain 



B 



\Wll 



dgpe 



J dt 

Foe d_ 

n dt 



B 



\Wll 



R 



B 



(eC-fieSB^^)) - 



m, 



Iwul^ilL — ■0) 



(110) 



This completes our solution for ge, as (110) is a closed equation for the part of gpe that 
is even in a. It is useful to note that only the zonal (i.e., a- independent) parts of C 
and 5B\\ will give a signiflcant drive for gpe - the a-dependent parts will tend to cancel 
under the flux-surface average. 
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4- 4- Field Equations 



In this section, we finally close our system of equations by writing the field equations in 
terms of our solution for 6fe and ion quantities. 

We find 6ip' by using the solution (81) for 6fe in the quasineutrality condition (29): 



rip ^ 



rflnA^e 

dip 



+ — [d^wge + Y.— fd^'wiK)^. (Ill) 



Similarly, we find 5B\\ by using (81) in (32) to find 



5B\\B 
- An 



^V±V±: ms{w±Wj_hs)^ + neV']_{C -6(f') 

s=i 



d\n.Ne dhiTp 
T. — h (Te - ev^o) 



dip 



dip 



;ii2) 



These two equations determine 6ip' and SB\\. 

It is instructive to note that (111) and (112) are not naturally solved in the moving 
coordinate system. It is because of this coupling to the fixed coordinates (and the ion 
physics) that we require the mapping between the coordinate systems (from the solution 
of (67)-(69)). 



5. Summary of the Reduced Electron Model 

The set of electron equations we have derived is complicated: the electron time-scale 
has been removed at the expense of introducing more dependent variables and more 
equations, some integro-differential. In these equations, we have the nine unknowns ip, 
a, 6u\\e, C; 9te, dps, hi, 6ip' , and We have provided a complete set of nine equations 
to determine these quantities. These are: 

• (67) and (68) to evolve ip and a which define the moving coordinate system that 
moves with velocity Wefi given by (79), 

• (87) to determine Su\ie, 

• and the coupled set (94), (103), (110), and (24) to evolve (, gte, gpe, and hi with 
the constituent relations (111) and (112) determining 6(p' and SB\\. 

Together, these equations comprise a rigorously-derived model for ion-scale turbulence, 
accurate up to corrections of order 6. The only assumption we have had to make, beyond 
our orderings, is that any stochastic component of the magnetic field (if present) satisfies 
5B,,och/B < e6. 



6. The Collisional Limit 



The equations for 6fe in Section 3.5 and Section 4.4 are derived under the assumption 
that z/ge ~ w ~ ^ll'^thr In many cases of interest, the collision frequency is in fact much 
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larger than this: u^e ^ oo- We can take this hmit as a subsidiary expansion, subsidiary 
to our expansion in the mass ratio. + 

Formally, we introduce the parameter u* = UeeqR/ Vth^, where qi^ip) is the safety 
factor (see Section 7.3 of [ ]) and R the cylindrical radial coordinate (the major radius 
of the tokamak). As our initial ordering was i/ge ~ ^, our ordering for u* is u* ~ S; for 
our coUisional ordering we let u* ^ 6 and expand in 




(113) 



The results of this expansion will be presented in the following sections. 



6.1. Maxwellian Electrons 

In our coUisional expanasion, C[ge] is the dominant term in (92). Thus, up to corrections 
that are small in u/uee 1, Qe is given by 

C[9e] = 0. (114) 

Hence, is a perturbed Maxwellian to leading order: 

where Suh and 6Tf, are arbitrary functions of r - 6nH will shortly be shown to be the 
homogenous (i.e., a function only of ip) part of the electron density perturbation and 
STe the perturbed electron temperature. 

Now, we know from Section 3.1 that ge is a solution of (49). Substituting this 
into (49), we see that dnn/Ne and 5Te/Te are functions of ip only (the derivatives 
in (49) being taken at constant Sf.)- Thus, in the this limit, collisions trap and de- 
trap electrons rapidly enough that there is no longer a distinction between passing and 
trapped electrons (a-dependence being the signature of kinetic trapped electrons) - all 
electrons sample the entire flux surface on the fluctuation timescale. 

Let us now further simplify our solution for g^. Substituting (115) into (81), we 
obtain 

^ t ) - - ^) ^ - (I - 2) - ° • 

From this expression, we observe that we were justified in our notation for 5Tf. - it really 
is the peturbed electron temperature relative to the moving surfaces. We also note that 
C and driH both contribute to the perturbed density. As we have the freedom to add an 
arbitrary function of ip to the definition of C - it does not change the velocity of the flux 

+ Strictly speaking, to treat this as a subsidiary expansion we also require that v^.^ <C fc||fthei so that 
it does not change the form of (43) and (46). However, if v^e is large enough that i^ee ^ ^||Wtiie, then 
the homogenous equation for is w\\b- Vge = C[5e]- The only solution of this equation is a perturbed 
Maxwellian with perturbed densities and temperatures that are functions of ip only. Thus, the results 
of this section still hold. 
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surface (see the end of Section 3.2), we redefine C to be C — TeSnn/eNe and so eliminate 
Shh entirely. Thus, our solution for 6fe becomes 

where the entire density perturbation is contained in the Boltzmann-like first term. 

6.2. Fluctuating Continuity Equation 

We now derive an equation for ( in the collisional limit. 

Using the solution (115) for ge in (94), we obtain, to lowest order in S/u*, 

d 



dt 



Up, + n — —- — < , -z— > rip \ i 



Te Tp T, B J ^ dl B eB y B ' Tp } T, 
dlnNp d\nTp\ \ dlxiNp 



+ l|{c>f}-^-vf-+(^^ + K.)-v(V-C)f 

, ^ xD^eA^lniVe d6B^inpTpd\nTp lupdud.- 



;ii8) 



da V eB J dijj da eB dip B dip dl 

where we have used the fact that, as we have absorbed the density perturbation into 
C, fd^wge = 0((5/z/*)ene). As in Section 4.1, this is just the fluctuating continuity 
equation for electrons. 

6.3. Fluctuating Energy Equation 

To close our system of collisional equations, we now need the fluctuating energy equation 
from which we will determine STp. Thus, we multiply (92) by meW^/(2Te) — 3/2 and 
integrate over all velocities to find 

dt 



" ~2 



35Te (r Tp6B\\ 6Tp 



l(7V. + 3V,).v(^-V^)^(f + lnT.)n. 
-Vn-vU-ip) ^^n, + Vo-V(V-C)^ 

^ ^ dip J-e 

dSBucTpUe fd In Np d\nTp\ lupdcud,- 



(119) 



da eB \ dip dip J B dip dl 

where q\\e = Jd^w {nipw'^ /2Te — 3/2) h'^\ As we do not know gye, we wish to eliminate 
it from this equation. Therefore, we average (119) over the perturbed flux surfaces and 



The Electron Response 



27 




T, 



■ V(V>-^)— ^ne) --(Vd- V(^-V')— y^^e) (120) 

where we have used (G.6) to ehminate some terms. 

Let us examine the terms in (120). Other than the last term on the right-hand 
side, which is due to the viscous heating of electrons, these terms seem fairly opaque. 
However, they all vanish if we apply the low-Mach- number ordering from Section 11.1 
of [i] - we can move rig outside the flux-surface averages and then use the properties 
of the fiux-surface average and the drift velocities to eliminate the terms in question. 
Thus, we conclude that these terms, which disappear if the Mach number is small, are 
to do with the exchange of potential energy of the electrons with thermal energy. 

Equation (120) completes our solution for 5fe- In the next section, we close the 
system by working out the electron contributions to the field equations in the collisional 
limit. 



6.4- Field Equations 

We now use our solution for 5fe to simplify the field equations. Substituting the solution 
(117) into (111) gives the following equation for 

S = l S = l 

Similarly, by substituting (117) into (112), we obtain an equation for 5B\\ in the 
collisional limit: 

5B B f 
-Vl^^ = ^^^^ • / rn,{w^w^h,)^ + rieVl (C - V) 

s=i 

+ (l-^)n.Vl^T.. 

The field equations (121) and (122) close the collisional electron model, which we 
now summarise. 



T 



dlnNe 

dip 



+ {Te 



c/lnTg 



(122) 



6.5. The Collisional Electron Model 

The collisional electron model is a simplified version of the model given in Section 5. 
We have replaced the distribution functions Qpe and gte by a single perturbed electron 
density and temperature. The complete set of equations describing this model is: 
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• (67) and (68) to evolve ip and a which define the moving coordinate system, 

• (87) to determine 6u\\e, 

• C is given by (118), 

• 6Te is given by (119), 

• (121) and (122) determine the fields 6ip' and 6B\\ respectively, 

• and the ion dynamics {hi) are given by (24). 

7. The Transport Timescale 

The procedure of Section 5 is sufficient to calculate the evolution of our plasma on the 
fast (turbulent) timescale. Our expansion in 6 also has consequences for the evolution of 
He, Te and u on the slow (transport) timescale. In Section 8 of [1], transport equations 
for He and were derived by considering the transport of particles and energy across 
ip = const, surfaces. It is clear from the results of Section 3.2 and Section 3.3 that, for 
electron transport, it is more appropriate to consider ip = const, surfaces. 

To this end, in Appendix H, we rewrite the electron kinetic equation relative to 
the moving ffux surfaces. We then average this equation over the turbulence and 
over the ffuctuating ffux surfaces. Finally, we integrate over velocity space to ffnd 
transport equations for particles and multiply by Eg and integrate over velocities to 
obtain the electron heat transport equation. This procedure is entirely analogous to 
that undertaken in Section 8 save for the extra complexity introduced by the rapid 
ffuctuation of the ip = const, surfaces. We summarize the results of these calculations 
in the following sections. 

7.1. Particle Transport 

In Appendix H.3, we derive the electron particle transport equation, correct to lowest 
order in 6: 



Note that whilst we use the same notation for the particle ffux Fe as in [ ], it is not the 
same quantity. In this paper the transport equations are derived for transport across a 
ip = const, surface not a ip = const, surface. 

Several important features of this ffux immediatly present themselves. Firstly, only 
the trapped electrons contribute - the passing particles play no part in the transport. 
Secondly, there is no "ffutter" transport (transport due to corellations between 6u and 
5-B|l despite the fact that the contribution to (170) from SA\\ is clearly not zero. This 
is, in fact, the primary advantage of writing our transport equations with respect to the 




(123) 



where S'e is deffned by (164) and the particle ffux is given by 




(124) 
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moving surfaces, we observe this cancellation analytically rather than having to recover 
it via a very careful evaluation of (166) and (170). Finally, we see that drifts that do 
cause transport are the drift due to ( and the V-B drift due to SB\\. This supports 
our interpretation of the drift duee to ( being the drift of electrons across the moving 
surfaces. 



7.1.1. Electron Transport in the Collisional Model As we proved in Section 115, in 
the collisional limit the electrons are trapped and detrapped fast enough to remove the 
distinction between trapped and passing particles - in effect all electrons are passing. 
This immediatly implies that if the electrons are collisional, there is no transport of 
electrons. Thus, the electron transport time is 0(5^^ (w/z/gg)) longer than the ion 
transport time - collisions are actually beneficial. Naturally, with sufficiently large 
collisions, neoclassical and classical electron transport reappear in the fluxes and the 
confinement worsens again. 



7.2. Electron Heat Transport 



Similarly, in Section Appendix 1.4, we derive the electron heat transport equation, to 
lowest order in 6: 

^5||^_rK),T.+ l|rte),= (125) 

pturb _|_ pcomp _|_ ppot _|_ (j{E) _|_ ^{E) 

where Ci^^ and Si^^ are defined by (192) and (193) respectively, the heat fiux is given 
by 



and the heat sources are given by 

p;omp ^ ^^^^ . y^^^ ^ 



turb/ 



(126) 



(127) 



ppot 



6(^0 



d_ 
dt 



eC 



turb / ^/) 



and 



P' 



turb 



d_ 
dt 



(eC-/ie5i?||) 



+ e(ne(CV-Wcff)turb) 



turb/ ?/) 



turb / t/i 



In these expressions for the heating, the divergence of u^s is given by 



V ■ u. 



off 



d_ 
dt 



5B\\ 

^l),aj B 



I{ip) duj c 
~l3~d^ ~ B 



(b • Vb) ■ 6 X V (C - V) • 



(128) 



(129) 



(130) 
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This expression is correct and closed to the required order in e and 5. Again, it is crucial 
to note that whilst we use the same notation for q^, P™™^, -Pf°*, and P^™^ as in [1], 
these are not the same quantities. 

Let us now interpret the heat flux and the heating terms. The heat flux parallels 
the particle flux: there is no flutter transport, only trapped particles contribute to the 
flux, and they contribute via the E x fi-drift across flux surfaces and the V-B-drift due 

to 

Of the heating terms, the flrst, P™™p is just the compressional heating due to the 
mean motion of the flux surfaces. The second, is the change in the thermal energy 
of the electrons due to the change in their electrostatic potential energy.* Finally, we 
have the turbulent heat source P*"''^. We interpret the terms in the expression (129) for 
P^™^ as follows: the first term is heating due to work done against the parallel electric 
field (through (); the term involving the time-derivative of SB\\ is in fact composed 
of the electron viscous heating, compressional heating due to the fiuctuating motion 
of the ip = const, surfaces and work done bending the magnetic field. This can be 
seen by substituting for 6B\\ from (130) - the term on the left-hand side of (130) is 
the compressional heating, the second term on the right is part of the electron viscous 
heating. The final term is the potential energy contribution to the compressional heating 
due to the fiuctuating motion of the exact fiux surfaces. 

7.2.1. Electron Heating in the Collisional Model Again, in the collisional limit, the fiux 
in (125) vanishes. We are just left with local heating terms, of which only the turbulent 
heating simplifies: 



P 



turb 




_ ) +e(ne(CV-iiefr)t,,b)v;' (131) 

/ turb/ ^ 

where V ■ itefr is given by (130) as before. Thus, despite the lack of any transport, the 
electron response to ion-scale turbulence can still dissipate free energy. All other terms 
in (125) remain unchanged. 

7.3. Momentum Transport 

As the momentum transport involves all species, not merely electrons, we do not write 
a new equation for it in the moving frame. The transport equation for toroidal angular 
momentum is (Equation (179) of [']) 

^||/'M« + ^|V-(.S?')^HS'"'>,. (132) 

where the moment of inertia of the plasma is now given by 

J = Y,^s{nsR^)^. (133) 



* For further discussion of the nature of this potential-energy-exchange term, see Section 8.3.5 of [1]. 
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and the total momentum flux is 

ySfX = {4**% + E + '"■-('/') («'r.>,) 



— dMihsW ■ V^)/A),„^, ■ R'V<j) (134) 



s=i 
1 



where the sums over s = i are taken over all ion species and fluxes vri'^"^-' and Tg for the 
ions are given by (184) and (170) of [1], respectively. The electron momentum flux is 
now (see Appendix H.5) 

(vr(^^))^ = (^m(^jd?w [m^wl - \m^wl^ g^l{^ - V^)^^ . (135) 

The three evolution equations (123), (125) and (132) close our system on the 
transport timescale. The fluxes are written entirely in terms of fluctuating quantities 
that our reduced electron model provides (note that nowhere do we require the part of 
gpe that is odd in a). 

7.3.1. The Collisional Limit In the collisional limit, Qe is Maxwellian, and so isotropic 
in velocity space - see (115). Therefore, the velocity integral in (135) vanishes and so 
there is no angular momentum transport due to the electrons. 



8. Summary 

In this section we summarise our work, flrst by comparing it in detail to previous 
models of the electron response, then by briefly collecting our results and outlining 
future applications of this work. 



8.1. Comparison to Previous Models 

Many previous works have introduced physically-motivated model responses for the 
electrons. Some are based on the linear electron response [ID] and others use fluid 
moments of the electron distribution function to try and capture the passing particle 
response [11, 14, 25]. Our model improves upon this by being a rigorous consequence of 
the electron gyrokinetic equation, clearly demonstrating how the electron response flts 
into the full turbulent transport framework. Also, by retaining the bounce- and flux- 
surface-averaged kinetics of the electrons, we retain the distinction between the passing 
and trapped populations and retain drift resonances that may be lost in a fluid picture. 
We recover some of the results of the fluid models in the collisional limit, where there 
are no trapped electrons. In particular, the fact that 6Th is only a function of ip means 
that the temperature has equillibrated along perturbed fleld lines as expected. Even in 
this limit, the previous models did not allow for a rapid toroidal rotation (comparable 
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to the ion thermal speed). Also all these models neglect 5B\\ and only retain 5A\\ in the 
perturbed magnetic field. 

One of the original and still widely-used models for the electron response is the 
adiabatic electron repsonse corrected for zonal perturbations [26] . It can be shown that 
this is a rigorous limit of our model in the limit of vanishing (3 (the electrostatic limit) 
and coUisional electrons. For coUisionless electrostatic turbulence, similar equations have 
been derived rigorously in [9]. Our model extends this to include fully electromagnetic 
turbulence and rapid toroidal rotation. The electrostatic model of [ ] also neither 
handles general magnetic geometry nor discusses the behaviour of passing electrons. 
This limit is useful as an easily-implemented stepping stone between the adiabatic 
electron model (stricly only valid for quite coUisional plasmas) and simulating the full 
electron gyrokinetic equation (which is either computationally prohibitive or formally 
invalid depending on whether one chooses to resolve the electron scales or not). 

The most complete previous model is that of [12]. However, this model is derived 
under two restrictive assumptions. Firstly, the mass ratio is assumed to be much smaller 
than we assume here (5 ~ e rather than 1 ^ 6 ^ e which we assume) and secondly 5B\\ 
is assumed to be identically zero. We note that (81) corresponds directly to equation 
(54) of [ ], but that the further assumption made in [12] of (in our notation) gpe = is 
unjustified {gpe = is not a solution of (110)) and hence neglects the electron response 
that is zonal with respect to the perturbed fiux surfaces. In addition, the model in [12] 
does not include the effects of sonic toroidal rotation. 

Finally, we must discuss how our model compares to the common practice 
of simulating the full electron gyrokinetic equation at ion scales (so-called "kinetic 
electrons"). Firstly, the electron gyrokinetic equation clearly contains electron-gyroscale 
and electron-bounce-frequency physics so simulating only the ion scales is formally 
invalid - deliberatly leaving the simulation unresolved in some sense. Secondly, how 
an unresolved simulation of the electron gyrokinetic equation behaves will be sensitive 
to the numerical algorithm used - it is possible to design a numerical algorithm that 
numerically bounce-averages the kinetic equation in the long-timestep limit, but of 
the commonly-used gyrokinetic codes only GS2 [27] does this. Indeed, the only way 
to confirm that simulating "kinetic electrons" gives the correct physics is to either 
compare with full two-scale simulations that resolve the electron dynamics as well 
as the ion dynamics or to compare with simulations of the electron model presented 
in this paper. As well-resolved two-scale simulations are generally beyond current 
computational resources, simulations of our model are the only way to confirm the 
validity of this practice. 

8.2. Conclusions 

In this paper we have derived equations describing the electron response to ion-scale 
turbulence. The fast electron timescales have been eliminated from this system. 
The scale separation between electron and ion scales both motivates and enables the 
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derivation of such a model. Starting from electron gyrokinetics ['], we systematically 
expand our equations in the square root of the electron-to-ion mass ratio to obtain our 
model. As discussed in the introduction, as studies of turbulence at non-zero (3 become 
more and more important - both for conventional and spherical tokamaks - having an 
electron response model that handles this rigorously will be necessary. We have derived 
full (/3 ~ 1) coUisional and arbitrary-coUisionality models - summarised in Section 6 
and Section 5 respectively. 

The procedure by which we constructed this model, a rigorous expansions in 
\fmjrni ^ 1, allowed us to clearly delineate which effects are due to electron-scale 
physics and which to ion-scale. Most importantly, the elimination of the electron scales 
has a profound impact upon the structure of the fluctuating magnetic field. We prove 
in Section 3.2, that purely ion-scale turbulence conserveres the total magnetic fiux, and, 
moreover that a velocity field can be found into which the magnetic field is frozen. 
This generalization of the usual MHD frozen-in theorem demonstrates that it is only 
electron-scale effects which can change the magnetic topology. The inclusion of /3 ~ 1, 
trapped particles, sonic rotation, and arbitrary axisymmetric geometry do not change 
this fundamental fact. 

Similarly, whilst deriving transport equations across the exact fiux surfaces, we 
have obtained another general property of ion-scale turbulence: passing electrons do 
not contribute to transport fiuxes (see Section 7). This result could prove useful in 
designing magnetic geometries that have beneficial trapped-particle properties, safe in 
knowledge that the passing particles never contribute to the transport. 

In addition to these general results, two other ancillary results have been proved in 
the course of this work. Firstly, in Appendix C, we proved that the topology of the mean 
magnetic field (i.e., the g(^/') profile) evolves on the resistive rather than the transport 
timescale. Importantly, this result holds independently of the nature of the turbulence. 
Secondly, in Appendix J, we prove that the commonly-used adiabatic electron model is, 
in fact, the correct model for coUisional [y ^ uS) and electrostatic (/3 <^ 1) electrons. 

Throughout our derivation, we have simply assumed that the magnetic stochasticity 
produced by electron-scale fiuctuations is sufficiently small that we can neglect it. This 
assumption must be investigated in detail and will be explored in future work. Similarly, 
the potential effects of small resonant regions around rational surfaces (which would 
give a finite width to each ■^/'-constant surface) should be examined closely. Indeed, a 
natural extension of this work would be to include the reaction of a small amoun of 
purely electron-scale turbulence upon the ion scales. As always, the final arbiter of 
the usefulness of these models will be their numerical implementation and application 
to simulating moderate- and high-/? ion-scale turbulence, in e.g. nonlinear studies of 
finite-/? ITG turbulence or kinetic ballooning mode turbulence. 
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Appendix A. Derivation of (42) 

The gyrokinetic potential x is 



x(r) = 6ip'(r) - —6Au(r) - -w± ■ 6A±(r). 

c c 



Using (16), we perform a Taylor expansion about r = Re to find 

^ M||(i?e) 



i^ll „ . , _ , w\\ b X w 



- -w±_ ■ 5A±{Re) - - 
c c 



b X w 



d 



dRe 

-SAJR, 



(Al) 



(A.2) 



fig <9-Re 

where we have used the fact that w ■ 6A/c ^ 6ip'. Gyroaveraging this expression, the 
third and fourth terms vanish to leave 



= {Re) - -^5A\\[Re) — + 0[5 x) 



(A.3) 



where we have used {w±wj_)j^ = w\/2{\ — bb), with I the unit dyadic, and 6B^\b 
V X (5A_L. 

Appendix B. Derivation of (57) 



Expressing E in terms of potentials, we arrive at the following expression for E ■ b: 



Eb 



IdA 



■b 



ld6A 



■b + 



c dt c dt 
Manipulating the first term of this expression, we find that 

b V^ = -uj(ij, t)b ■Vij = --b-(Bxu), 

c c 

where we have used (6), (9), and (12). Using B = B — SB, we obtain 



5. V$ = 



1 



B - 6Bj X u 
u ■ (V(5A) b-b- (V(5A) • u 
u-V (6A-Pj -b-V{5A-u 
u ■ y5A\ 



6 ■ [(V X 5 A) X u] 



b-V {6A-u) + -u- (V(5A) ■ 6B 
B 



uVb-bVu 

+ (e 



5A 



B 



(B.l) 
(B.2) 



(B.3) 
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where in the last hne we have used (45), u ■ Vb = b ■ Vu ((A. 17) of [I]), and neglected 
6B ■ [Vu) ■ 6A as it is 0(e^ (vi-^./c) B). Substituting (B.3) back into (B.l), we obtain 



cB _ 

The final term in this expression is 0(e'^ {vthjc) B) and can be neglected. We can then 
use 6(p' = 6ip — 6 A ■ u/c to finally write the parallel electric field as 

Eb=-bV^-b- W --^■b--(^ + u-v] 6A:\ + O ( e^'^B) . (B.5) 

c ot c \ot J V c / 

Appendix C. Derivation of (58) and Conservation of Mean Magnetic Flux 

In this Appendix, we first derive (58), as required for the derivation in Section 3.2. 
Then we prove that there is a field-line-presrving velocity for the mean magnetic field, 
irrespective of the fiuctuations. Finally, we demonstrate that this implies that there are 
two timescales for the evolution of the mean magnetic field: the transport timescale te 
on which the fiux ^ipi^R, z, t) evolves and the resistive timescale r,, ~ 5~^te on which the 
safety factor q{ip,t) evolves. 

We now proceed to derive (58) in a way that exactly parallels the derivation of (54) 
from (24). The neoclassical drift-kinetic equation that determines Fie is ((112) of [ ]) 

^W^-gi;- (C[Fie])^ = -Voe ■ ^ + ^^ll^oe^ ■ b + (C[FoJ)^. (C.l) 

The first term on the left-hand side and the second term on the right-hand side are 
0(e^5~^f2jFoe), all other terms are O(e^fijFoe). Thus, to lowest order in 6, we obtain 

6 dA 

. VFi, = —W^lFoe^ . b, (C.2) 

where we have dropped the distinction between and r as it is small in e. This 
equation is completely similar to the first line of (46) and so we seek solutions in an 
identical fashion. Let (c.f. (48)) 

Fie = A(r, Ee, He, Cr)Foe + Ge{r, Be, /ig) + O (c^Foe) , (C.3) 

where Ge satisfies 

w\\b-VGe = Q. (C.4) 

As Ge is axisymmetric and so cannot depend on a, we immediatly solve (C.4) to find 
Ge = Ge{ip,ee, fie)- Substituting (C.3) back into (C.2) and dividing by Foe, we obtain 

p 3 A 

u,6.VA = 6. (C.5) 
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By the same logic as in Section 3.1, A is a function of r only (this time there is no gradient 
term to drive a temperature perturbation). Moreover, by absorbing any density moment 
of Ge into A, we can assume that jd^wGe = and so 

A = ^ (C.6) 

where nie is the first-order correction to the mean electron density: riie = jd^wFie. 
Using this result in (C.5), we finally obtain 

5.v!^ = ff.6. (C.7) 

He TeC at 

This is precisely (58) as required in Section 3.2. 

We now go to next order in 6 to obtain an equation for Gg- To this order, we have 

w\\V\\F['J - C[Ge] = -VDe ■ VV'^, (C.8) 
where F^l^ is the next correction to Fig. Using (116) of [ ], we have 

Multiplying this equation by Ge/FQ^, integrating over all velocities, and flux-surface 
averaging, we obtain 

d^w^G[Ge]) =0. (C.IO) 

The only solution to this equation is for Gf. to be a Maxwellian with a density and 
temperature that are functions of t/^ only. We can absorb any such solution into Foe- 
Thus, Ge vanishes and the solution for Fig is merely 

i^ie = — i^Oe + 0(e(5Fg) (C.ll) 

Substituting (C.7) into the expression for E ■ B vn terms of the potentials and 
A, we have 

^ . B ^ . (v. + if ) ^ -B . V (. + ^) + O {f,^-fB^) . (C.12) 

This form of the mean electric field is entirely analogous to (55). Thus, by the arguments 
of the second paragraph of Section 3.2, the mean magnetic flux is conserved and there 
is a field-line-preserving velocity for the mean magnetic field. 

Finally, we look at the implications of this for the evolution of the mean magnetic 
field given by (6). This field is determined by the two functions /('?/',t) and ■?/'(/?, 2;, t), 
which a priori we expect to evolve on the transport timescale te- In [ ], it is shown 
that instead of /, we can use the safety factor q{;i\}) given by 

q{^) = l^y'm{R-'),. (C.13) 
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to characterise the magnetic field. It is this quantity that will evolve slowly in the limit 
of 5 ^ 1. The evolution of q{ip) is given by [ ]: 



di 



Using the expression (C.12) for E ■ B, we see that 

{E.B).^ = 0[e's'^B'), (C.15) 

and so the right hand side of (C.14) is 0{e'^6Qi) = 0(r~^), and thus q evolves on 
the resistive timescale, not the transport timescale. Since q{ip) is the number of times 
a field line winds around the vertical symmetry axis for one poloidal transit around 
the magnetic axis, a topological quantity, it changes only when the magnetic topology 
changes. However, the shape of the fiux surfaces changes on the transport timescale 
in response to changes in the pressure profile (via the Grad-Shafranov equation, see 
Section 7.3). It is therefore clear from (C.13) that /(V") also changes on the transport 
timescale. 

Appendix D. Representations of 6B 

Once we introduce the Clebsch form of B we have two equivalent expressions for 6B: 
firstly, the usual expression in terms of 6A\\ and 

6B = bx V±M|| + 6Biib (D.l) 

and secondly the expression in terms of and a, 

6B = ^-4) xVa -Vtfj X Va. (D.2) 

Let us investigte this second expression for 6B. Writing the gradients of and a in 
terms of the (■?/',«,/) coordinates we have 

Examining the size of the various terms in this equation we see that 

^ - 1 ~ 0(1) and ^ - 1 ~ 0(1), (D.4) 

because ip and a can vary on the short perpendicular spatial scale (V_l ~ Pi'^)- Thus, 
the first term in (D.3) is comparable to B. Similarly, we have that 

9 'i,-A^0(.t), ,D.5) 



dl dl \ ) \ (i- , 
as parallel variation is much weaker than perpendicular. Thus, the second and third 
terms in (D.3) are smaller than the first by one power of e. In order that bB is indeed 
smaller than B we require that 
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Thus, to lowest order, ip and a are not "independent". This represents the fact that, to 
lowest order, ■UeS is an incompressible flow and thus cannot generate any SB^\. Taking 
the scalar product of (D.2) with Vip we obtain 

84' 



6B-V^ = -B 



dl 



Substituting for 5B via (D.l), we are able to relate ip to 5A\\: 

d5A\\ _ di) 
da dl 

Similarly, we can find an expression for 5A\\ in terms of a: 
d5A\\ _ _da 
dip ~ ~~dl' 



(D.7) 



(D.8) 



(D.9) 



Hence, given any one of 6A\\ ,ip or a we can determine the other two from the appropriate 
pair of equations from (D.6), (D.8) and (D.9). We do not use (D.3) to determine 6B\\ 
since for this we must know tp and a to O(e^) or equivalently the O(e^Vthi) compressible 
corrections to Weff ^ and we do not. Instead we use (32) to determine 6B\\ since it 
involves only lowest order quantities. 

From the above results we can also form two expressions for V(5y4|| in terms of a 
and ip: 



V±M|, = - — {a - a)Vij + —{i! - ij)Va, 



(D.IO) 



and 



d ^ d ~ 

-—Ja - a)Vi) + -~(^ - i))Va, 
dl dl 



VxM|| = -^~{a - a)Wip + ^~{ip - ip)Va, (D.ll) 
We will require these expressions when we come to eliminate 6A^\ from our equations. 



Appendix E. Derivation of (79) 

In the definition (56) of we expand E in terms of potentials, keeping terms up to 
0{evthi), to obtain 



UcS =— I 1 



SBii 



b X V ($ + <^o + (^V^ - 



{u ■ V6A) X b 



B - ' (E.l) 

+ f/(r)b + 0(eV)- 

where we have used the fact that (9/(9t + u ■ V) (5A ~ efij(5A to rewrite the time 
derivative of 8 A in the second term. Using (12) in the first term in (E.l), we find 



B 



6B\\ 
If 



6 X V$ 



6B\\ 
If 
6B\\ 
~lf 



u 



u 



B 
B 



+ ^SB^ X ViP 



B 
6B± 

~l3~ 



(E.2) 



wb + 0(eW 
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where we have also used (6) and (45). The second tern in (E.l) can be rewritten as 
follows: 

4 (w • VSA) X 6 = 4 [(V(5A) ■ w - w X (V X SA)] x b 



1 ~ 1 ~ 

--b X V (5A ■ w) - - (w X X b + 0(e2t;th, 
B B 



(E.3) 



-—b X V {5A ■ u) 
B 



SB, 



b ■ u 



u 



B B 

Substituting (E.2) and (E.3) back into (E.l), we obtain 

Ues =u+^bxV (<^o + V - + \u{r) -b-u\b + 0{e\i,J, (E.4) 
B L J 

where we have used b ■ u = lui^p, t)/B and (45). Letting U{r) = b ■ u and substituting 
for ^ from (60), we find 



u^s = u + —bxV [ 6lp' 



IniV, 



erif. 



e +0(eV), (E.5) 



where we have used the fact that {c/B)b x V (Tenie/ene) ~ e^'^thi and can thus be 
neglected. Finally, picking ^ = lnA'^e('^); "we arrive at 

c 



u^s = u + —bxV 



6(p' H iip — ip 



e \ / dtp eUe 



+ 0(eV)- 



(E.6) 



This is precisely (79), in which ( is defined by (80). 



Appendix F. Properties of the Bounce Average 

In this appendix we prove the properties of the bounce average that were stated in 
Section 4.2. 

First, as the second portion of the integration only differs from the first in the sign 
of a and the direction of integration, we have that, for any function a that is independent 
of 0", 



w\\a{l,tlj,a,ee,iJ,e,t) 



0. 



(F.l) 



Sceondly, combining (65) and (99), we discover that, for any single-valued function a, 

da' 



w\\b ■ Va(/, ijj, a, Se, jJ^e, o", t) 



w 



dl 



dl 



\W\\ 



II 

da{a 



dl 



da{a = — r 
dl 



(F.2) 



0, 



where the contributions from the endpoints has vanished as at the bounce points 
a(cr = +1) = a{(x = —1) in order that a is continuous through = 0. 

We now prove that (Vue ■ VV^^ = 0. Starting from the formula for V^e, (89), we 
can show that, to leading order in 6, 

Vne = ^ V|,^,^^,, X , (F.3) 
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where the gradient is taken at constant energy Se and magnetic moment /Xe- Using this 
form for Vde, we have that 

VDe-V7/i= ^V- (u;||bx V^) . (F.4) 
By using the expression of a divergence of a vector field A: 



V-A = 

we can write (F.4) as 
d 




w^\B 
fip 1 da 



B 



6 X W 



) ■ Va 


d 




+ —~ 






dl 


B 



where we have used (F.2). Now, 

i (6 X Vip) ■Va = bb = l + 0(e^ 
B ^ ' 



Thus, 



m(,cw\\ d 



e dliB 



W\\ 



B V 



VB 



(F.5) 



(F.6) 



(F.7) 



(F.8) 



as w\\ is axisymmetric and so independent of a and thus a (to lowest order). We can 
now take the bounce average of this expression and use (F.2) to get the required result. 

Finally, we prove that the bounce average commutes with the time derivative in 
the field-aligned coordinates to the required order in e. The integration over / is taken 
at constant t, ip, a, Eg and fig so clearly commutes with the time derivative in these 
variables. The fact that the denominator in (99) varies slowly in time, follows from the 
fact that it is in fact just the time taken for an electron with a given location, energy Se 
and magnetic moment fie to travel along a field line from one bounce-point, to the other, 
and back again. As the 6Bj_ ~ eB and ~ /, the contribution to the distance along 
the perturbed field line from the fluctuations is small. As w\\ varies slowly in space, its 
value along the perturbed field line is almost the same as its value along the unperturbed 
field line. Thus, the bounce time along the perturbed field is a slowly varying function. 



Appendix G. Properties of the Flux-Surface Average 

As the flux surface average is defined in exactly the same way as the average over 
equilibrium flux surfaces defined in [1] (the fact that the integral is taken at constant Se 
and He rather than v does not affect the derivation), we can use the results of Section 3.3 
of [1] with a simple change of notation {ip becomes ip). Thus, for any function A{r,v) 
we have 

{Ah = - !^A=^ f^^A, (G.l) 
where the integration is taken over the surface labelled by ip and 

V = hm i / d'r = [ (G.2) 

Ai,^0 Alp J A J B 
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(V-A)^ = -— (V'(A-VV^) J, (G.3) 



Similarly, for any vector field A(r, v) we have 

v di) V ' ^ 

where the derivative are all taken at constant Ee-, /ie and d. 

We now prove the results that are stated without proof in Section 4.3. Using (G.3), 
we can show that 

'b . Va)^ = (V . {ab) )^ = 1 1 [V {aB . V^)^) = 0. (G.4) 

Using (G.l), we see that for any single valued function a{r) 

to lowest order in e, where we have used B = B + 0{eB), axisymmetry, and periodicity 
in a. For any fluctuating function a{r) and any fluctuating function we have that, 
to lowest order in e, 

5 , ,A / B^ da\ db Id fB^a\\ db ,^ , 

where we have used (G.5) and (66). 

We now prove some results on fiux-surface averages involving Voe and Vd- First, 
multiplying (F.8) by B/\w\\ \ and using (G.4) gives 

/AvDe-V7/;\ =0. (G.7) 

\klll 

Secondly, we prove (108). Working from (F.3), we have 



\w 



^VDe-VA\ =-l—V-{\w\\\bxV\) 



(G.8) 



eme V dip 

However, as, to lowest order in e, B and are independent of a, we have 

\w\\ \bxVX- Vip)_ = {\w^\ \BX) ^ _ , (G.9) 

where we have used the fact that A is a fluctuating quantity to neglect dX/dl compared 
to dX/da. Inserting this into (G.8) and using (G.5), we obtain (108) as promised. 
Finally, as we can write Vd in the form 

^^ = — Vxf-lV (G.IO) 



e \B 

we can prove that 



Vd-ViP). = 0, (G.ll) 



The Electron Response 



42 



and 



Vd- VA 



0, 



(G.12) 



in the same manner as (G.7) and (108), by using the fact that Tg = Te('?/') + O(eTe) 
which commutes with the flux-surface average. 

As the magnetic perturbation is small, the fluctuations cannot change the area of 
a flux surface by more than an a fracation of order 0(e). Thus V can change by at 
most eV in one fluctuation turnover time, and so dV /dt ^ uV . This means that the 
flux-surface average commutes with the time derivative at constant if), <5, /: 



dt 



(G.13) 



Appendix H. Transport Equations for Electrons 

In this Appendix we derive the transport equations presented without proof in Section 7. 
We proceed from the kinetic equation written in r and v variables 



d \ e 

— +U-V \ fe + 'W-Vfe+ — 

ot I nip 



E + -V X B 

c 



dfe 

dw 



It will be convenient to introduce the pecualiar velocity w 
we can write (H.l) as 



C[f4+Se. (H.l) 

V — Wcff; at which point 



dfe 

dt 



+ V|~-(We)+lt^-V|~/e + 

eB 



d_ 
dw 

dfe 



fe 



e du^s ^ 

V ■ VUeff 

TTLf, dt 



(H.2) 



H w X b 

crrie dw 



C[fe]+Se 



To explicitly evaluate fluxes, we will need to use (H.2) to evaluate the gyrophase- 
dependent piece of fe correct to second order in e. Taking (H.2), dropping terms smaller 
than 0(e^r2i/e), and dividing by eB/cnie, we have 

dfe 



b X w ■ 



cm. 



dw eB 



W-Vfe+{ ~ Ve 



cnip 



-V ■ \/u 



+ 



cnie 
eB 



d_ 
dt 



B 



+ Woff ■ V 



cB 

fe-C[fe 



cflf 



dfe 

dw 



(H.3) 



In this equation, the second-order parts of fe appear only on the left-hand side and so 
we will use this to eliminate the second-order parts of fe from our calculation. 

The structure of this appendix is as follows: flrst, in Appendix H.l and 
Appendix H.2 we introduce some averages that will be needed later. In Appendix H.3, 
we derive the particle transport equation for electrons across the moving surfaces. In 
Appendix H.4, we derive the transport equation for electron heat. Finally, we evaluate 
the flux of toroidal angular momentum due to the electrons. 

Detailed derivations of the intermediate results used in this appendix are provided 
in Appendix I. 
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Appendix H.l. The Turbulence Average at constant ip: 

In order to derive the transport equations, we will need some new averages. We define 
a new average over the fluctuations where the perpendicular spatial integrals are taken 
between surfaces of constant ip and a rather than ip and a. Similarly, the time integral 
in this average is taken at constant ijj,a and /. We denote this new average by (■)tui.b- 

r/2 

iaYturk = j dij j da j dtg / T j dij j da, (H.4) 

T/2 ' 

where the spatial integral is taken over a small patch that has dimension A (the 
intermediate length scale) in each direction - for reference compare with (15) of [l]. 

This average commutes with time derivatives at constant ip and with flux-surface 
averages over -^-constant surfaces. 

As usual, the composition of this average and the average over ^/'-constant surfaces 
is an average over time composed with an average over an annular region between ip = i/jq 
and = V'o + It is important to note that this average differs from the turbulence 
average in [ .] only by terms that are small in e. This will be crucial, as we wish to 
replace all instances of this average by their unprimed version in the final equations. 

Appendix H.2. Flux Surface Averaging on the Transport Timescale: 

We will want to average the transport equations over the fluctuating flux surfaces. For 
this we need two identies. The first is just a specific case of (G.3). If A(r, t) is a function 
only of space and time then 

where the derivative on the right is taken at constant t only. 

Secondly we need to interchange flux-surface averages and time derivatives. We do 
this via 

'dg\ _ I d , , Id 

m/^-y'dt 

which is proven in an identical fashion to (45) of ['] except that we retain the velocity 
explicitly rather than replacing it with dip/dt via (67). 



^ V dip \ ' i> 



Appendix H.3. Particle Transport 

We now move to the meat of the problem. Deriving our first transport equation. 
Integrating (H.2) over all velocties 

d_ 



j d'wf, + V ■ d^wf.w^ + (V ■ u,s) J d'wfe = Si^\ (H.7) 
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where we have used the definition of the time derivative at constant ip, a and /. 
Averaging this equation over the fiuctuating fiux surfaces we obtain 



1. ^ 
V' di 



V' 



(H.8) 



where we have used (H.5) and (H.6). We now muhiply this equation by V', apply the 
new average over the fiuctuations, and divide by V to find 

V' di 



(H.9) 



where we have used the fact that iV) = V'{tp) = V + 0{eV') and the fact that 

\ / turb 

:(") 



Ss ~ rie/TE- The particle fiux Fe is defined by 



(re), 



I turb 



(H.IO) 



In Appendix 1.3, we use the kinetic equation to find the following explicit expression 
for the particle fiux: 



turb/ ip 



(H.ii; 



Appendix H.4- Energy Transport and Heating 



Multiplying (H.2) by msw'^/2, integrating over all velocities and averaging over fiux- 
surfaces and the turbulence we obtain 
1 3 d 
'V'2 di 



cff 



dt 



turb 



where the electron thermal energy fiux is defined by 



d'^w-rriew'^ feW ■ \/tp 



^/ turb 

We show in Appendix 1.4 that we can rewrite this equation as 
1 3 d 
V'2 di ^ 

pturb _|_ pcomp _|_ ppot _|_ (j{E) _j_ 

The heating terms are, in order, the turbulent heating 

d 



^ V'dip 



p: 



turb 



d W{ Qe 



dt 



(eC - /Xe5fi||) 



turb / {j) 



the mean compressional heating 

pcomp ^ ^^^^ . y^^_ 



(H.12) 



(H.13) 



(H.14) 



+ e(ne(CV-lieff)turbh, (H.15) 



(H.16) 
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where is the velocity of the mean flux surfaces defined in (44) of [1], and the exchange 
between thermal energy and electrostatic potential energy: 



d_ 
dt 



In (H.17) and (H.15), the divergence of u^s is given by 



(H.17) 



V ■ Ucff = 777 
Ot 



i.iT + (S ■ V^) ~ - I (f ■ V6) . 5 X V (C - V) . (H.18) 



This expression is correct and closed to the required order in e and 5. 

In Appendix 1.5, we show that the heat flux can be written in terms of ge as 

Mi = (/ (^."e^'^ - e^„ - ec) (9,4 (< + i|^B,) , (H.19) 

Appendix H. 5. The Electron Angular Momentum Flux 

We start from (C.15) of ['] and drop all terms that are small in 5 to obtain 

(V^) ■ n, ■ (i?2V0) = /m,i?" I (fw{{w ■ V0) [{w X 5B) ■ V0] 5f,)\ , (H.20) 

\ J I turb 

where we have expanded ba^ and only kept terms that are of order e^R^BuiTi, i.e., 
comparable to the ion momentum flux. Substituting for Sfg from (81) and performing 
the gyroaverage gives (134). 

Appendix I. Derivations for Appendix H 

Appendix I.l. The solution for /g 

We now collate results from the rest of the paper to write down a convenient form for 
/e. First, using the general solution (13) for /g, we can write 

/e =Foe(V^(i?;),£e) f 1 + + /^^(^^(ii'J, «(i?'J, /i„ ^, t) 



J-e J (1.1) 
+ Fie{R'e, ^e, /^e) + /2e, 

where /2e ~ e^/e and we have defined h'^ = he — [ipi^R'^) — ip{R'e)]dF()e/dip in order to 
be able to evaluate Fq^ at ipiR'e)- In addition, we have introduced the guiding centre 
position defined with respect to the exact field 

R'=r- crrie 1.2 

cB ^ ' 

which differs from R'^ by 0{e'^R'J. 
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We now use the solutions for the distribution function obtained in this paper. For 
Fie, we have from (C.ll), that 



rile 



Foe + F^^\r,ee,fie) + 0{6\''f\ 



(1.3) 



where F^^ ~ 6efe contains the higher-order (in 6) contributions to Fie. We also wish 
to write in terms of the new velocity variable w. To this end, we introduce the 
following quantity 



F' 

e 



rrie 



3/2 



exp 



27rTe(^(i?'J) 

me (UeS -U) -W 



+ 



meUj\ij{K))R^ 



rripW 



Te(^(il'J) 2Te(V^(i?'J) 2re(V^(i?'J) 



(1.4) 



+ 



Te 



Foe- 



Using (60) and (80), we obtain the following expression for ^ 
— = InNeitp) + ^ + — + 7f^, 

J- p J- p. ^p J- p 



(1.5) 



in which Tg is evaluated at and we have let ^ = 0. Substituting this into (1.4), and 
using the definition (21) of Foe, we can show that 



Tie Te 



+ FI 



(H) 
e 



( lJ^' ( H) 

where F{^ ~ 6efe differs from by some small gyrophase-independent pieces - we 
will never need F{^ explicitly so can absorb any gyrophase-independent function into 
it. Because of the convenience of the form of F^, we choose to write fe as 



, e (C - Sip') ,j^{H)',,,f 

Je - ^e ^ ^Oe + ^ le + + he- 



(1.7) 



Using our solution (81) for 6fe and the definition of h'^, we have that, to lowest order in 



n-e - ^Oe + Qe- 

^ e 

We will use this equation later to eliminate h'^ from our final results. 



Appendix T2. An equation for f2e 

As mentioned above, we will need properties of f2e to calculate the fluxes in our transport 
equations. We now derive such an equation by using the solution (1.7) in (H.3) - we 
will need to retain corrections to f2e up to 0(5e^/e). We perform this substitution step- 
by-step. The preliminary step is to notice that all the terms in the second line of (H.3) 
are 0(5^e^/e) or smaller and so do not contrbiute in this calculation. 
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Next, we compute the left-hand-side of (H.3) for F^: 

of: 



w X b 



dw 



crrie _ ^~ \ dTe 
— —w ■ Vip < — - 

cB [d^ _ 



where we have used the definition of and the fact that h ■ Vip = 0. Using this result 
and the definition (1.4) of F^, we can show that 

(wxb)-^-^w-\/F^ = -^w-(eV^+lm,u\^)VRAF^ 

CTfi Tfl 
+ -^W • V (C - V) + -T^WW : V (iteff - It) Fq^, 

eB -Le 

where we have used (79) to show that {ucs — u)-bx w = {c/ B)w-V (C — ^v')- However, 
we also have that 

c „^ crrie „ \ dF'„ cnie ^ ( „^ 1 9, ,.„„9\ ^ 

- ■ Viteff • = - ■ eVe + -m,u\ij)VR^ F 

cnie rrie ^ ^ — ^, 

^— : Vu^sF^ 

eB J^e 

where we have used the fact that v = u^s + w and the fact that, to lowest order in e, 
Ucs = u. Thus, the contributions from F'^ to (H.3) cancel (up to corrections of order 
0(5^e^/e)) save for the last term in (1.11). 

Then, using these results and substituting (1.7) into (H.3), we have 

(r ~\ df2e Cme ( ^ ^ \ C , % 

I' " "^J ■ 5^ =^ r ■ - ■ ^0 + + ^) ■ 9^ ^2) 

eB J-e 

where we have used (1.8) to write the result in terms of Qe and used the fact that, for 
any function H{R'^,ee, fie), 

bxw) ■ — = ^wVl H. 1.13 
/ dw eB 

We now convert the derivative of Qe at constant w into a derivative at constant Se 
and fie- 

-W-Vge= ^W- VL^.„^ ge ^W ■ £e77- + VL fi 



e 

e 



eB " eB eB \ dSe 

cme ^1 crue f ^ % 95(e\ 



eB ^^'^^ eB \ de^ dfie 



(1.15) 



e 
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in which we have used the definition Eg = {me/2)w'^ — cipo + 0(5^Te). Combining this 
equation with the term from (1.12) involving ipo, we have 

CTTle^ C dge 

eB B dw 

crrie ~ ^ cme / _ _ , e _ — \ dg, 

By taking the derivative of the definition of /ie we can show that 

w ■ V Us — ew± ■ Vipo = —rrieW ■ In B + -^'"^ " — ew^ ■ Vipo (1-16) 

We can rewrite this in terms of Vbe as 
w ■ Vfj,s — ew± ■ Vipo = 

TJleW^- i^bxVBe) ^ iw ± i V b + ■ V \n B ^ 

Substituting this back into (1.15) and thence back into (1.12) we arrive at: 

bxw)- -— = ^VC • 77^ ^—ww : (Vu) Foe 

J dw B dw eB Te 



crrip 



nieW ■ & X Voe ^ ( w^wx_ i Vb + -^b ■ V In 5 j 



dge (1.18) 



cB 

Estimating f2e from this equation, we see that f2e ~ ^e^/e, rather than f2e ~ e^/e which 
we would naively expect - this is because it is only the gyrophase-dependent piece of 
/2e we need and, to lowest order in 5, f2e is gyrophase-independent. 

Appendix 1.3. Derivation of (H.ll) 

Using the solution (1.7) for /g in the definition of the particle fiux, we have 

V'Ve = l^j d'w - ^^^^^Foe + K + f2e) W ■ V^^ ) . (1.19) 

\ ^ ^ I turb 

Substituting for from (1.4) and performing the velocity integration, we have 



V'Te=(v'(^ne (iteff - *t) ■ VV^)^ 



+ d'wiK + f2e)w-ViP 



turb , , 

(1.20) 



turb 

Using (79) for Ugff and using the definition of the fiux-surface average in the first term, 
we have 



turb 

(1.21) 



^'r.^(/...&|(C-V) 

+ y' i^j d''w{h'^ + f2e)w-V^^ 



turb 
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Integrating by parts with respect to a in the first term, gives 



turb 

(1.22) 



'^1 turb 

The calculation of the term involving h!^ is particularly involved. In Appendix 1.6, 
we prove the following identity: 



/ dadl^iC - Sip')— " 
J B oo! \ 



B 



(1.23) 



turb 



Using this in (1.22), we find 



"I" I turb 



^ \' ^^^^ X "^i^ ■ w (w X b 



dw 



(1.24) 



V*/ turb 

where we have integrated by parts in the term involving f2e- 

Now, the derivative of f2e in (1-24) is precisely the one appearing in (1.18). Thus, 
substituting (1.18) into (1.24), we have 

d f w'\ 

''VI turb 

+ {V' { I d?w^ (bxVC- V^a) I (1.25) 

/ turb 

where we have integrated by parts with respect to w in the second line and used the 
gyrophase- independence of in the third line. We have also used the isotropy of Fq^ 
and the gyrotropy of ge to eliminate terms in (1.18) that contain an even power of 'W± 
(we are multiplying by one power of w± so this becomes an odd power and thus vanishes 
upon integration). 
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We now prove that the last hnes of (1.25) are small and can be dropped. Writing 
it in terms of /ig and using (F.8), we obtain 

■;2 



C 



Wll 



'I'/ turb 



B 



dge 



turb 



(1.26) 



v 



d^wwnb -Vl 



turb 



0, 



where we have been able to write (F.8) in terms of w because the difference between w 
and w is 0(efthj and thus negligible. 

Finally, integrating by parts with respect to a in the second line of (1.25) and 
dropping all the references to lu, we have 



(1.27) 

turb / 'tp 

where the passing electron response has vanished because dgpe/da = (at constant 
and fie)- This is (H.ll) as required. 



Appendix 1.4- Derivation of the Heat Transport Equation 

As in Appendix C.4 of [1], it is convenient to adjust the form of the pressure 
evolution equation before calculating the heat flux. In particular, we subtract 
(V ■ [e{ipo + C) / d?'^'^fe\)^ from both sides of (H.12), and use 



turb 



turb 



j d^WeV ifo ■ (iy/e)turb + ^'^oV ■ d^Wwfe 



turb 



d^weVipo ■ wfe ) 

I turb 



(1.28) 



- ev^o 



d_ 
Ft 



\ J / turb 
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where we have used (H.7). This results in the following heat transport equation 



1 3 d 
^2 di 



1 



du^s 



turb / ip 



eV (e - ^o) 



dt 



V ■ VWcff 



^1 



turb 



^/ turb 



(1.29) 



+ ( ev^o 



d_ 
dt 



_ne + /(V- Wcff) jd^wf. 



5'(") 



turb 



Appendix I.^-l- Electron Heating. We now evaluate the heat sources in the temperature 
evolution equation, from left to right. The first heating term is 



i(^(^jd^wf,w-V{^-y^o] 



turb 



d^wfeW ■ VC 



, (1.30) 



turb 



where we have used (1.5) for C,- Rearranging the derivative to act on /e, we have 
' d^wfeW ■ VC 



turb 



e( ( V ■ d^wCwf, 



(1.31) 



turb 



/ turb 

Now, integrating (H.3) multiplied by {eB/cme)C over w, we have that 

d 



= C / d-'wW ■ V| ~ /e + C^eV ■ Weff + C / d'w 



dt 



fe 



(1.32) 



where we have integrated by parts with respect to w as needed. Substituting this back 
into (1.31), we have 



d^wfeW ■ VC 



turb 



e ( V- d'wCwf, 



turb 



A - 



d 



+ ^<<^^^^-"^^Wturb-^\\^^9t 



(1.33) 



c 



V"' turb 

where we have integrated by parts with respect to time in the second term. The first 
of these terms cancels with the corresponding term in (1.29) and the second and third 
will become part of the turbulent heating. We will need to process terms involving 
V ■ Weff several times in this section. In Appendix 1.4.2, we show that we can write it in 
a closed (if unwieldly) form - see (1.45). For conciseness we will not expand V ■ itefr in 
this section. 
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Now, the time-derivative term is small: 



rfWe^-ti}/e) ) =0{6e^QiT,), (1.34) 

and can thus be neglected. 

Substituting the solution (1.7) for /e into the remaining heating term gives 

- mJl^jd'wv ■ Viieff • w (f^ + K + ^-^^^Foe + f[^^' + h)j y) . (1.35) 

\ ^ ^' I turb 

Estimating the size of the contribution from the gyrophase-dependent piece of 
/2e ~ f^e^/e, we see that it is 0(5e^fijneTe) and so we can neglect it. Now, as 
(V ■ itofr)turb ~ ■ Vitcff • b\ ~ O(e^fij), we can also neglect the contributions from 





turb 
{ H)' 

and any gyrophase-independent piece of f2e- Of the remaining terms, we handle 
the Fg term first. Using (1.4), we have 

—niei { / d^WV ■ VWefl ■ wF' 



H turb (1-36) 

= - ((V ■ U^s)turh^e)^ Te - me (((Ueff " VUeff) ■ (WcflF - w))turb^e)^ , 

where we have used the fact that, to lowest order in e, the flux-surface average and the 
turbulence average commute. As ■Ucfr — it ~ ei^th^ we can use (79) to express u^s in 
the second term. Doing this, we find that the entire second term is 0{6e^neTeQi) or 
smaller and can thus be neglected. In the first term, we cannot use (79), but instead, 
we interchange the divergence and the turbulence average to find that 



-mJ{ / d^WV ■ VWeff ■ wF^ ) ) = -Te (rieV ■ V^)r , (1.37) 




where we have used {Ues)^^^^.^^ = - the average of the velocity of the exact surfaces 
is just the velocity of the average surfaces. We now work out the corresponding term 
involving Qf,. Using the gyrophase- independence of Qe, we have 




d wv ■ Vwcff ■ 'i^g, 

/ turb 

/ d'^wB (^1 - bb^ : (VWeff) f^e9, 

' V-/ turb 

Substituting for V ■ ■Uefr from (1.45), we obtain 

turb 




(1.38) 



me( ( jd wv- Vweff • wge ) ) = - { { Id Wfieg^ 





(1.39) 



V"' turb 
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Finally, we have to deal with the term involving cpo and the divergence of u^s- 
Using our solution for /g, we have 



e ( (V ■ u, 



eS 



turb / t/i 



e {(foUeV ■ V^)^ + e{(po{ { — rie + jd wg^ ] V ■ Ues 



turb/ 



Using (1.45) to substitute for V ■ u^s in the last term, we have 



e ( Vol (V ■ u. 



turb / t/i 



+e ( (/9o< 



d_ 
dt 



5B\\ 



bb : VwcflF 



turb / ip 



where the parallel compression is given by (1.46). 

Using (1.33), (1.37), (1.39), and (1.41) in (1.29), we obtain 



turb / i/j 



d^wge 



d_ 
dt 



(eC + fXeSBn) 



d_ 
dt 



He + UeV ■ + 



•>p,ct,l 



^' turb 

^ne+ [d:'wg)jV-u,s) - Si^^ 

J J I turb 



(1.40) 



(1.41) 



(1.42) 



-T,(neV-V^)^ + Cf) + 5P, 

where V ■ ■Uefr is a shorthand for the expression given in (1.45). Gathering the terms 
in this equation together according to their physical interpretation, we obtain the heat 
transport equation (H.14) as required. 



Appendix 1.4-2. The Divergence of Ue^: We wish to calculate the divergence of ■UeS 
correct to 0{e'^Qi). We start from Faraday's Law: 

— = V X (u,s xB)+ Oie^SBQi). 

Taking the component of this equation along b, we obtain 
d 



(1.43) 



+ Weff ■ V 5Bn =B{\-bb] : Vueff + 0(e'5Bni). 



dt 

Rearranging this into an equation for V ■ u^s, we have 



V ■ u 



eff 



Ft 



5B\ 



+ bb : Vu,s + 0(e^6BQi). 



(1.44) 



(1.45) 



ip,a,l B 

To the required order, we can use (79) to write the parallel compression due to u^s as 
bb : Vu,s = (b ■ V^) - I (b ■ Vb) ■ b X V (C - V) + 0{e^6Q,). (1.46) 
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Appendix 1.5. Derivation of (H.19) 

The derivation of the heat flux proceeds in exactly the same way as the particle flux 
derivation. 

Using the solution (1.7) for /e in the definition of the heat flux, we have 
V'q,= lv'(^j d^w{^m,w^-e^o-edj {F, + h'^ + h,)w -V^"^ \ . (1.47) 

\ ^ I turb 

Substituting for F'^ and integrating, we have 



turb 



' "^1 turb 



(1.48) 



where, for conciseness, we have defined Be = {l/2)meW^ — e(po. 

As with the calculation of the particle flux, we use the identity (1.23) to handle the 
term involving h'^. Using (1.23) in (1.48) and integrating by parts where required, we 
obtain 



turb / 



+ (v' Q d^weji xVi'-wlwxby^'^ (1.49) 
- (v'leCfeW ■ Vij' 



'^1 turb 

where we have used (79) and the definition of the flux-surface average to cancel the first 
term in (1.48) with the contribution arising from the first term of (1.23). 
Substituting for /2e from (1.18), we have 

turb 



turb 



^{^'{I '^^'^^'^ -Vcibxw- V^) Qe^ _ ^ (1.50) 
-(v'l^Jd'w'j^^i^V^^.Vi^g^ 



turb 



5 



turb 

where we have integrated by parts with respect to w and used the gyrophase- 
independence of ge- We have also used the isotropy of Fq^, and the gyrotropy of g^. 
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to show that most terms in (1.18) do not contribute to the heat flux. The proof that 
the term involving Vde can be dropped is identical to the proof for the particle flux - 
the manipulations of (1.26) go through in exactly the same way with an extra multiple 
of Ee- 

Next, we show that the third and fifth lines of (1.50) cancel. Taking the fifth line 
and expressing /e via (1.7), we have 

V'(^j d^weCfeW ■ VV^^ ^ = (y'i^j d^'^eC (F^ + h'J w ■ V^pj ^ , (1.51) 

where we have been able to drop terms associated with F^^ as it is both small and 
gyrophase- independent. Substituting for from (1.4), we have 



V'i / d^weC 



TJlf. (■Ucflf ~'")''"^n I b X w 



w ■ Vip 

■>pl turb 



where we have expanded h'^ around r = R'^ and dropped terms that are small in S. 
Substituting for h'^ from (1.8) and using (79), we see that the first term in the brackets 
in (1.52) cancels with the contribution from the first term on the right-hand side of (1.8) 
leaving us with 



turb 



V'{ I «eC( -Vgew-V^ 

i'^ turb 

Performing the gyroaverage explicitly we have, 

7,\ \ _ / t7'/ f ^3~._ /-..a^Se 



(1.53) 



da 



turb \ ^ / turb 



(1.54) 



turb 

where we have used the definition of the flux-surface average to integrate by parts with 
respect to a. Now, taking the third line of (1.50) and performing the gyroaverage, we 
obtain 



turb 



^/ turb 

Thus, the third and fifth lines of (1.50) cancel. 



, , (1-55) 
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Finally, using this result and integrating by parts with respect to a in the second 
line of (1.50) and dropping all the references to lu, we have 

qe=(f d'wUe (Im^w' - cifo) ^ U + |f ^^nl ) ) , (1-56) 



where the passing electron response has vanished because = (at constant Eg and 
fie)- This is (H.19) as required. 



turb/ 1/1 

— n fci 

da 



Appendix 1.6. Derivation of (1.23) 

In order to derive this result, we will need a concise notation for the following operation: 

V^ = V-V/~— - = Va— + V^— (1.57) 

dl da ^ dilj 

This is an important operation as V/ is not parallel to the exact field, but the strong 
perpendicular variation is contained in and a. 
Now, the quantity we have to evaluate is 

(v'(K,^.Vi)}^^^^j^d((^^.^.i)^^ (1,58) 
First, we evaluate B{ip,a,l) at ip{R'g), a{R'^), and / to obtain 



/ turb 



(1,59) 



In the term involving the gradient of In B we can perform the gyroaverage to find 
j dadl(^L{h'^w±)^-Vip^ 

\ \ ^ 1 r I turb \ / turb 

We see that this term is just the V-B drift in the exact magnetic field. 

Now we handle the first term on the right-hand side of (1.60). Integrating by parts 
with respect to space, we have: 



turb 

(1.60) 



/dadll 

^ ,r ,turb \ ^ 



dadl{ ( -J^ — ) ■ VV' ) = / dadll w± -V 



— / dadll ( w I -V 



(I.6I; 



e 



Sc^ie B(Ri., 

r I turb 



where we have used the fact that Vip = V±ip and the fact that we can then choose to 
take that derivative at constant Ee and fie rather than at constant w. 
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The crucial step in this derivation is to now write inside the first bracket as 
ipiR'e) + {r — i?g) ■ V±ijj. Doing this, we obtain 

hi ^ 



dadli 



B(R' 



turb 



dadl(^± -Vj 



h' 



£e,Me \B(R' 



~ , b X w 



(1.62) 



turb / 



— / dadli ( w±_ V _ 



h' 



ee,l^e B(R' 



4j 



turb 



Now, for any function H{ip{R'^),a{R'^),l,ee, fie) we have that 

dRl dH\ ^ /diV 



wx_ -Vj 



H 



dR' 



0. 



(1.63) 



Applying this result to (1.62), we obtain 



dadli 



BiRL 



J dadl(^ 



turb 



h' 



b X w 



(1.64) 



\B(R' 



turb/ 



We now work on writing this term as a perpendicular divergence. First, we convert the 
derivative back to one at constant w: 



dadli { -^J^ 
\ \ B{R'^ 

V' i B( wVi 



K 



B(RL' 



turb 

b X w 



(1.65) 



turb I rl ij) 



-V 



2VL, 



b X W 



turb/ ?/) 



Next, using the fact that d{\^^^/da ~ ea and {dSe/da) ~ B{dfie/da) ~ eTg by 
axisymmettry, we have: 

hi 



dadli 



BiRL 



V' ( BV 



B(R' 



turb 

b X w 

~n7~ 



(1.66) 



turb I r I ij) 

We can now perform the gyroaverage inside the divergence to find that 

h' 



dadli 



B(R' 



-V' (^BV_ 



r 



turb 



^' " ' b X 



fi2a 



/ turb / r / ^/j 



(1.67) 



The Electron Response 58 
j dadlUj^w\ -vA = 0((5eV/e|VV5|i?). (1.68) 



BiR'e) / / u 

^ 1^' I r I turb 

Inserting this final result back into (1.25), we obtain 

j dad!(^l{Kw^)^-V^^ = j d5^dll^^^h'}j . (1.69) 

Substituting for h'^ from (1.8) and using the fact that dB/da = d5B\\/da. we obtain 
(1.23) as required. 

Appendix J. The Electrostatic Limit 

In this Appendix we provide the electrostatic limit of our equations. 
Appendix J.l. The Electrostatic Ordering 

Formally, this is an expansion in <^ 1 subsidiary to our expansions in 5 and e. 
In the main body of the paper the Alfven velocity f a is assumed to be comparable to 
the ion thermal velocity fthj we must now distinguish the two. We assume that the 
fluctuations we are interested in have frequencies comparable to those of sound waves, 
rather than Alfven waves: 
d 



and that typical fluctuating velocities scale with the ion thermal speed (rather than the 
Alfven speed): 

j d^whsW ~ evtY,,ns, (J. 2) 

and that this scaling also applies to the fluctuating E x B velocity: 

Q 

— bxV6ip^evth,. (J.3) 
Jd 

In order to consider electrostatic turbulence, we assume that 6B± and 6A\i are 
small: 

JD C 

These orderings imply that ip and a differ from ip and a by amounts that are small in 

Estimating the size of terms in (112) we see that 

^ ~ (J.5) 

and so we can drop all terms containing 6B\^ . 

As we wish to obtain the simplest equations possible, we will also neglect all finite- 
Mach-number effects. Thus, we combine this expansion with the low-Mach-number 
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expansion of Section 11 of [^]. The primary result of the low-Mach expansion is that 
we can neglect ipo and the poloidal variation of the mean density. In addition to this 
simplification, we can drop the distinction between S(p' and 6(f - from now on only 
working with the electrostatic potential dip. 

In the subsequent sections we will apply this ordering to our coUisionless and 
coUisional equations. One result is common to both the coUisional and coUisionless 
electrostatic limits - the general solution for (. We obtain this by applying our 
electrostatic ordering to (82). The left-hand-side of (82) is ^/P smaller than the right; 
thus, to lowest order, we have 

b ■ V (C - 5<^) = 0, (J.6) 

where we have also used the fact that b = b + 0(a//3) Solving this for ( we have 

c = 5^+m, (J.7) 

with ( an arbitrary function to be determined. Using this result in the expression for 
the field-line velocity Ueff , we see that Ues is directed within the flux surfaces and merely 
moves one field line onto another without distorting them. Thus, this is consistent with 
our electrostatic assumption that 6B is small. 

Appendix J. 2. The CoUisional Electrostatic Limit 

In the coUisional limit, we have already used the freedom in the definition of ( to 
eliminate 6nH and so cannot use this freedom to eliminate (. Instead, we use (118) to 
determine (. Dropping all terms that are small in ^/^, (118) becomes: 

d 



dt 



{S^}^ + C ) = 0, (J.9) 



where we have used the fact that ip — ip e^/]3^|J. Using (J.7) for C and averaging over 
a flux surface, we obtain 

d_ 
dt 

where we have used (G.ll) and (G.12). Assuming no constant perturbations, we have 

C = -{5^)^. (J.io) 

Applying the low-Mach- number limit of [ ] to (120), we obtain 
'5Th /25Bi\\ \ 2 / 1 d ,~ ,A.rfw 



where we have used the low-Mach-number result that Ue is a flux function to remove it 
from the flux-surface averages whereupon we used (G.ll) and (G.12) to eliminate most 
of the terms in (120). Now, using (J. 5) and (J. 4) in (J.ll) and again assuming that 
there are no constant perturbations, we have that 6Th = 0. 
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Substituting these results into (117) we have 



Sfe = ^ ^-Foe + O(ev^Foe), (J.12) 



with 6ip found from 




e 



U^-{5v)\ + ^Y.zJSw{K)^. (J.13) 

s=i 



This is just the usual adiabatic electron model including the zero response to zonal 
perturbations [ ]. This reassures us that we can recover known physics from our 
complex model and that there is a physically-relevant regime where the adiabatic 
electron model is an exact limit of gyrokinetics and not an ad-hoc prescription. 

Appendix J. 3. The Weakly- Collisional Electrostatic Limit 

In this limit, we can use the freedom in the definition of C to set = in (J. 7). Using 
(79) we see that this implies that 

Wcff = w + 0(et;th,y/3). (J. 14) 

Equation (103) for gte then becomes 

+ ^ ■ gte + (Vce ■ Vgte)\\ + ^ - {{C[gte])R)\\ 

+ U-W \ ^ Foe ' - 




dt J Te da dip ' 

where we have used (J.4) and (J.5) to drop terms involving 6A\\ and SB\\. This is the 
usual bounce-averaged kinetic equation originally derived in [9], again demonstrating 
that the preexisting exact models are contained in our model. However, we must still 
find a solution for the distribution function of passing electrons: gpe (which was not 
discussed in ["]). This is determined from the electrostatic limit of (110): 

From this equation we see that gpe is driven by the competition between the shielding 
effects of the passing electrons that drive gpe towards the maximally-shielding solution 

1 /BS^\ eF„,_ 

and collisions which drive gpe back towards a Maxwellian. The solution for gpe is non-zero 
and so the original discussion in [ ] is incomplete. Thus the solution for 5fe is 

5fe=^-^+ge, (J.18) 
with ge = gpe + gte found from (J. 15) and (J. 16) and Sip obtained from 
1 / Z'^n^e e \ ^ 1 /" .■j 




^e \ Ta Te 



6ip = - [d'wge + y— [d^w{hs),. (J.19) 

ne J Ue J 

s=i 
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